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ABSTRACT 

We report a measurement of the Type la supernova (SN la) rate in galaxy clusters at 0.9 < z < 1.46 from 
the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN la rate measurement 
with detected z > 0.9 SNe. Finding 8 ± 1 cluster SNe la, we determine a SN la rate of 0.50to.i9 ( stat ) -om 
(sys) hj SNuB (SNuB = 1(T 12 SNe L~ X B yr _1 ). In units of stellar mass, this translates to 0.361^3 (stat) 

io'.oe ( s y s ) SNuM (SNuM = 10" 12 SNe M~ 1 yr _1 ). This represents a factor of w 5 ± 2 increase over 
measurements of the cluster rate at z < 0.2. We parameterize the late-time SN la delay time distribution with 
a power law: Sk(t) oc t s . Under the approximation of a single-burst cluster formation redshift of Zf = 3, our 
rate measurement in combination with lower-redshift cluster SN la rates constrains s = — 1.41^o.4o» consistent 
with measurements of the delay time distribution in the field. This measurement is generally consistent with 
expectations for the "double degenerate" scenario and inconsistent with some models for the "single degen- 
erate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental 
dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red- 
sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. 
Finally, the upper limit of one host-less cluster SN la detected in the survey implies that the fraction of stars in 
the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts. 
Subject headings: Supernovae: general — white dwarfs — cosmology: observations 
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1. INTRODUCTION 

Type la supernovae (SNe la) are widely accepted to be 
the result of the thermonuclear explosion of a carbon-oxygen 
(CO) white dwarf (WD). The explosion is believed to occur 
as the WD nears the Chandrasekhar mass by accreting mass 
from its companion star in a binary system. Despite the con- 
fidence in this basic model, many uncertainties remain about 
the process that leads to SNe la (see Livio 2001, for a re- 
view). Chief amongst them is the nature of the companion 
donor star. The leading models fall into two classes: the single 
degenerate scenario (SD; Whelan & Iben 1973), and the dou- 
ble degenerate scenario (DD; Iben & Tutukov 1984; Webbink 
1984). In the SD scenario the companion is a red giant or 
main sequence star that overflows its Roche lobe. In the DD 
scenario, the companion is a second WD which merges with 
the primary after orbital decay due to the emission of gravita- 
tional radiation. 

A better understanding of the SN la progenitor is de- 
manded from both an astrophysical and a cosmological per- 
spective. Astrophysically, SNe la dominate the produc- 
tion of iron (e.g., Matteucci & Greggio 1986; Tsujimoto et al. 
1995; Thielemann et al. 1996) and provide energy feedback 
(Scannapieco et al. 2006) in galaxies. Knowledge of the 
SN la rate is necessary to include these effects in galaxy 
evolution models. However, an accurate prediction of the 
SN la rate in galaxies of varying ages, masses and star for- 
mation histories requires a good understanding of the nature 
of the progenitor. This is particularly true for higher red- 
shifts where direct SN rate constraints are unavailable. From 
a cosmological perspective, the progenitor has become a cen- 
tral concern following the use of SNe la as standardizable 
candles in the discovery of dark energy (Riess et al. 1998; 
Perlmutter et al. 1999). With hundreds of SNe now being 
used in the precision measurement of cosmological param- 
eters (e.g., Hicken et al. 2009; Amanullah et al. 2010), astro- 
physical sources of systematic error will soon become signif- 
icant. While the unknown nature of the SN progenitor system 
is unlikely to bias measurements at the current level of uncer- 
tainty (Yungelson & Livio 2000; Sarkar et al. 2008), it could 
become a significant source of uncertainty in the future, as it 
leaves open the question of whether high-redshift SNe are dif- 
ferent than low-redshift SNe in a way that affects the inferred 
distance. 

Measuring the SN la rate as a function of environ- 
ment has long been recognized as one of the few avail- 
able methods for probing the SN la progenitor (e.g., 
Ruiz-Lapuente et al. 1995; Ruiz-Lapuente & Canal 1998; 
Yungelson & Livio 2000). SN la rates constrain the progeni- 
tor scenario via the delay time distribution (DTD), where "de- 
lay time" refers to the time between star formation and SN la 
explosion. The DTD is the distribution of these times for a 
population of stars, and is equivalent to the SN la rate as a 
function of time after a burst of star formation. The delay 
time is governed by different physical mechanisms in the dif- 
ferent progenitor scenarios. For example, in the SD scenario, 
when the donor is a red giant star the delay time is set by the 
time the companion takes to evolve off the main sequence. In 
the DD scenario, it is dominated by the time the orbit takes 
to decay due to gravitational radiation. The result is that the 
shape of the DTD depends on the progenitor scenario. 

However, the interpretation of the DTD is complicated by 
its dependence on other factors, not all of which are com- 
pletely understood. These include the initial mass func- 



tion (IMF) of the stellar population, the distribution of ini- 
tial separation and mass ratio in binary systems, and the 
evolution of the binary through one or more common enve- 
lope (CE; see, e.g., Yungelson 2005) phases. Theoretical de- 
lay time distributions were computed analytically following 
the proposal of both the SD (Greggio & Renzini 1983) and 
DD (Tornambe & Matteucci 1986; Tornambe 1989) scenar- 
ios. Later, theoretical DTDs were extended to include var- 
ious subclasses of each model and a wider range of param- 
eters (Tutukov & Yungelson 1994; Yungelson & Livio 2000; 
Matteucci & Recchi 2001; Belczynski et al. 2005; Greggio 
2005). In various recent numerical simulations, different 
plausible prescriptions for the initial conditions and for the bi- 
nary evolution have lead to widely ranging DTDs, even within 
one scenario (Hachisu et al. 2008; Kobayashi & Nomoto 
2009; Ruiter et al. 2009; Mennekens et al. 2010). A measure- 
ment of the DTD then must constrain not only the relative 
contribution of various progenitor scenarios, but also the ini- 
tial conditions and CE phase, which is particularly poorly con- 
strained. Still, most simulations show a difference in the DTD 
shape between the SD and DD scenarios. In both scenarios, 
the SN rate is greatest shortly after star formation and gradu- 
ally decreases with time. However, the SD scenario typically 
shows a strong drop off in the SN rate at large delay times not 
seen in the DD scenario (but see Hachisu et al. 2008). 

The DTD can be measured empirically from the SN la rate 
in stellar populations of different ages. Measurements corre- 
lating SN rate with host star formation rate or star formation 
history have now confirmed that the delay time spans a wide 
range, from less than 100 Myr (e.g., Aubourg et al. 2008) to 
many Gyr (e.g., Schawinski 2009). Correlations with star 
formation rates (Mannucci et al. 2005, 2006; Sullivan et al. 
2006; Pritchet et al. 2008) show that SNe with progenitor ages 
< a few hundred Myr comprise perhaps ~50% of all SNe la. 
Measurements as a function of stellar age (Totani et al. 2008; 
Brandt et al. 2010), show that the rate declines with delay 
time as expected. 

It is more straightforward to extract the DTD in stellar pop- 
ulations with a narrow range of ages (with a single burst of star 
formation being the ideal). Galaxy clusters, which are dom- 
inated by early-type galaxies, provide an ideal environment 
for constraining the shape of the DTD at large delay times. 
Early-type galaxies are generally expected to have formed 
early (z > 2) with little star formation since (Stanford et al. 
1998; van Dokkum et al. 2001). Cluster early-type galaxies 
in particular form even earlier than those in the field, with 
most star formation occurring at z > 3 (Thomas et al. 2005; 
Sanchez-Blazquez et al. 2006; Gobat et al. 2008). Measuring 
the cluster SN la rate over a range of redshifts from z = 
to z > 1 provides a measurement of the SN la rate at delay 
times from ~2 to 11 Gyr. Obtaining an accurate rate at the 
highest-possible redshift is crucial for constraining the shape 
of the late-time DTD: a larger redshift range corresponds to a 
larger lever arm in delay time. 

In addition to DTD constraints, there are also strong moti- 
vations for measuring the cluster SN la rate from a perspec- 
tive of cluster studies. SNe la are an important source of iron 
in the intracluster medium (e.g., Loewenstein 2006). Cluster 
SN rates constrain the iron contribution from SNe and, paired 
with measured iron abundances, can also constrain possible 
enrichment mechanisms (Maoz & Gal- Yam 2004). The high- 
redshift cluster rate is particularly important: measurements 
show that most of the intracluster iron was produced at high 
redshift (Calura et al. 2007). The poorly-constrained high- 
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redshift cluster rate is one of the largest sources of uncer- 
tainty in constraining the metal-loss fraction from galaxies 
(Sivanandam et al. 2009). 

Cluster SNe la can also be used to trace the diffuse in- 
tracluster stellar component. Intracluster stars, bound to 
the cluster potential rather than individual galaxies, have 
been found to account for anywhere from 5% to 50% 
of the stellar mass in clusters (e.g., Ferguson et al. 1998; 
Feldmeier et al. 1998; Gonzalez et al. 2000; Feldmeier et al. 
2004; Lin&Mohr 2004; Zibetti et al. 2005; Gonzalez et al. 
2005; Kricketal. 2006; Mihos et al. 2005). The use of 
SNe la as tracers of this component was first demonstrated by 
Gal- Yam et al. (2003) who found two likely host-less SNe la 
out of a total of seven cluster SNe la in 0.06 < z < 0.19 
Abell clusters. After correcting for the greater detection ef- 
ficiency of host-less SNe, they determined that on average, 
the intracluster medium contained 20^2% °f me total cluster 
stellar mass. The intrinsic faintness of the light from intra- 
cluster stars, combined with (1 + z) A surface brightness dim- 
ming, makes surface brightness measurements impossible at 
redshifts much higher than z = 0.3. Type la supernovae, 
which are detectable up to and beyond z = 1, provide a way 
to measure the intracluster stellar component and its possible 
evolution with redshift. 

The cluster SN la rate has recently been measured at lower 
redshifts (z > 0.3) in several studies (Sharon et al. 2007; 
Mannucci et al. 2008; Dilday et al. 2010), and at intermedi- 
ate redshift (z ~ 0.6) by Sharon et al. (2010). However, at 
higher redshifts (z > 0.8), only weak constraints on the high- 
redshift cluster la rate exist, based on 1-2 SNe la at z = 0.83 
(Gal- Yam et al. 2002). In this paper, we calculate the SN la 
rate in 0.9 < z < 1.46 clusters observed in the HST Clus- 
ter Supernova Survey. We address the host-less SN la frac- 
tion, and use our result to place constraints on the late-time 
DTD in clusters. Maoz et al. (2010, hereafter MaozlO) have 
already combined our results with iron abundance measure- 
ments and rate measurements in other environments to place 
even tighter constraints on the SN la DTD. 

This paper is organized as follows. In §2 we review the sur- 
vey, placing particular emphasis on the aspects relevant to the 
rate calculation. In §3 we describe the selection of supernova 
candidates used in this rate calculation and the determination 
of supernova type for these candidates. In §4 we carry out 
efficiency studies to determine the detection efficiency of our 
SN selection. In §5 we measure the luminosity of the clusters 
based on data from the survey. In §6 we present results and 
characterize systematic errors. We discuss interpretations for 
the delay time distribution and conclude in §7. Throughout 
the paper we use a cosmology with Ho = 70 km s _1 Mpc -1 , 
£Im = 0.3, ft \ = 0.7. Unless otherwise noted, magnitudes 
are in the Vega system. 

This paper is one of a series of ten papers that report 
supernova results from the HST Cluster Supernova Survey 
(PI: Perlmutter, HST program GO- 10496), a survey to dis- 
cover and follow SNe la in very distant clusters. Paper I 
(Dawson et al. 2009, hereafter Dawson09) describes the sur- 
vey strategy and discoveries. This work, Paper II, reports 
on the SN la rate in clusters. Paper III (Meyers et al. 2011, 
hereafter Meyers 11) addresses the properties of the galax- 
ies that host SNe la. Paper IV (Ripoche et al. 2011) in- 
troduces a new technique to calibrate the zeropoint of the 
NICMOS camera at low counts rates, critical for placing 
NICMOS-observed SNe la on the Hubble diagram. Pa- 



per V (Suzuki etal. 2011) reports the SNe la lightcurves 
and cosmology from the HST Cluster SN Survey program. 
Paper VI (Barbary et al. 2011) reports on the volumetric 
field SN la rate. Melbourne et al. (2007), one of several 
unnumbered papers in the series, present a Keck adap- 
tive optics observation of a z = 1.31 SN la in i^-band. 
Barbary et al. (2009) report the discovery of the extraordi- 
nary luminous supernova, SN SCP06F6. Morokuma et al. 

(2010) presents the spectroscopic follow-up observations for 
SN candidates. Finally, Hsiao et al. (in preparation) de- 
velop techniques to remove problematic artifacts remaining 
after the standard STScI pipeline. A separate series of pa- 
pers, ten to date, reports on cluster studies from the survey: 
Hilton et al. (2007); Eisenhardt et al. (2008); Jee et al. (2009); 
Hilton et al. (2009); Huang et al. (2009); Rosati et al. (2009); 
Santos etal. (2009); Strazzullo et al. (2010); Brodwinetal. 

(2011) ; Jee et al. (2011). 

2. THE SURVEY 

The details of the HST Cluster SN Survey are described 
in Dawson09. Here, we briefly summarize the survey and 
highlight the details relevant to the rate calculation. The sur- 
vey targeted 25 massive galaxy clusters in a rolling search 
between July 2005 and December 2006. Clusters were se- 
lected from X-ray, optical and IR surveys and cover the red- 
shift range 0.9 < z < 1.46. Twenty-four of the clusters have 
spectroscopically confirmed redshifts and the remaining clus- 
ter has a photometric redshift estimate. Cluster positions, red- 
shifts and discovery methods are listed in Table 1 . Note that 
cluster positions differ slightly from those reported in Daw- 
son09 due to the use of an updated algorithm for determining 
cluster centers. 

During the survey, each cluster was observed once every 20 
to 26 days during its HST visibility window (typically four 
to seven months). Figure 1 shows the dates of visits to each 
cluster. Each visit consisted of four exposures in the F850LP 
filter (hereafter zs^o)- Most visits also included a fifth expo- 
sure in the F775W filter (hereafter 2775). We revisited clusters 
D, N, P, Q, R and Z towards the end of the survey when they 
became visible again. 

Immediately following each visit, the four z^5o expo- 
sures were cosmic ray-rejected and combined using Mul- 
tiDrizzle (Fruchter & Hook 2002; Koekemoer et al. 2002) 
and searched for supernovae. Following the technique em- 
ployed in the earliest Supernova Cosmology Project searches 
(Perlmutter et al. 1995, 1997), we used the initial visit as a 
reference image, flagged candidates with software and then 
considered them by eye. Likely supernovae were followed up 
spectroscopically using pre-scheduled time on the Keck, and 
Subaru telescopes and target-of-opportunity observations on 
VLT. For nearly all SN candidates, either a live SN spectrum 
or host galaxy spectrum was obtained. In many cases, spec- 
troscopy of cluster galaxies was obtained contemporaneously 
using slit masks. Candidates deemed likely to be at higher 
redshift (z > 1) were also observed with the NICMOS cam- 
era on HST, but these data are not used in this work. 

A number of visits were contingent on the existence of an 
active SN. At the end of a cluster's visibility window, the last 
two scheduled visits were cancelled if there was no live SN 
previously discovered. This is because a SN discovered on 
the rise in either of the last two visits could not be followed 
long enough to obtain a cosmologically useful light curve. In 
addition, supplementary visits between pre-scheduled visits 
were occasionally added to provide more complete light curve 
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Table 1 

Cluster positions and redshifts 



ID 


Cluster 


Redshift 


R.A. (J2000) 


Decl. (J2000) 


Discovery 


A 


XMMXCS T2215 9-1738 


1.46 


22 h 1 5 m 59 s 


— 17° 37' 59" 


X ray 


B 


XMMU J2205.8-0159 


1.12 


22 h 05 m 50 s . 6 


—01° 59' 30" 


X-ray 


c 


XMMU J 1229. 4+0151 


0.97 


12 h 29 m 29 s . 2 


+01° 51/ 21" 


X-ray 


D 


RCS J022 1.6-0347 


1.02 


02 h 21 m 42 s . 2 


—03° 21' 52" 


Optical 


E 


WARP J1415. 1+3612 


1.03 


14 h 15 m 11 s . 1 


+36° 1 2' 03" 




p 


TSCS T1432 4+3332 


1.11 


1 4 h S2 m 28 s 1 


[ uu •J'J \J\J 


TR-Snit7pr 


Q 


TSCS T1429 3+3437 


1.26 


14 h 29 m 17 s . 7 


IQ40 07/ 1 0// 


TR-Snit7pr 


H 


ISCS J1434.4+3426 


1.24 


14 h 34 m 28 s . 6 


+34° 26' 22" 


IR-Spitzer 


I 


ISCS J1432. 6+3436 


1.34 


14 h 32 m 38 s . 8 


+34° 36' 36" 


IR-Spitzer 


j 


TSCS T1434 7+3519 


1.37 


1 4 h 04m a as n 

It Ot: t:U .U 


+35° 1 Q' 36" 


TR-Snit7er 


K 


TSCS T1438 1+3414 


1.41 


14 h S8 m 08 s 2 

It UO • Zj 


+34° 14/ 13" 


TR-Snitzpr 


L 


TSCS T1433 8+3325 


1.37 


14 h 33 m 51 s . 1 


+33° 25' 50" 


TR-Snit7er 


M 


CI T 1604+4304 

V_x 1 J 1 \J\J \ \ \ *j\J'-x 


0.90 


1 fi h 04 m 2S S 8 


+4S° 04/ 37" 




N 


RCS T0220 9-0333 


1.03 


02 h 20 m 55 s 5 


—03° 33 ; 10" 




P 


RCS J0337.8-2844 


l.l a 


03 h 37 m 51 s .2 


-28° 44' 58" 


Optical 


Q 


RCS J0439.6-2904 


0.95 


04 h 39 m 37 s . 6 


-29° 05' 01" 


Optical 


R 


XLSS J0223.0-0436 


1.22 


02 h 23 m 03 s . 4 


-04° 36' 14" 


X-ray 


S 


RCS J2156.7-0448 


1.07 


21 h 56 m 42 s . 2 


-04° 48' 04" 


Optical 


T 


RCS J151 1.0+0903 


0.97 


15 h ll m 03 s .5 


+09° 03' 09" 


Optical 


U 


RCS J2345.4-3632 


1.04 


23 h 45 m 27 s . 2 


-36° 32' 49" 


Optical 


V 


RCS J2319.8+0038 


0.90 


23 h 19 m 53 s .4 


+00° 38' 13" 


Optical 


w 


RX J0848.9+4452 


1.26 


08 h 48 m 56 s . 4 


+44° 52' 00" 


X-ray 


X 


RDCS J0910+5422 


1.10 


gh 1Q m 45 s A 


+54° 22' 07" 


X-ray 


Y 


RDCS J1252.9-2927 


1.24 


12 h 52 m 54 s .4 


-29° 27' 17" 


X-ray 


Z 


XMMU J2235.3-2557 


1.39 


22 h 35 m 20 s . 8 


-25° 57' 39" 


X-ray 



References. — A (Stanford et al. 2006; Hilton et al. 2007); B,C (Bohringer et al. 2005; Santos et al. 2009); D (also known as RzCS 052; Andreon et al. 2008a,b); D, N, U (Gilbank 
et al. in prep); E (Perlman et al. 2002); F (Elston et al. 2006); G, I, J, L (Eisenhardt et al. 2008); L (Brodwin et al. in prep; Stanford et al. in prep); H (Brodwin et al. 2006); K 
(Stanford et al. 2005); M (Postman et al. 2001); Q (Cain et al. 2008); R (Andreon et al. 2005; Bremer et al. 2006); S (Hicks et al. 2008); V (Gilbank et al. 2008); W (Rosati et al. 1999); 
X (Stanford et al. 2002); Y (Rosati et al. 2004); Z (Mullis et al. 2005; Rosati et al. 2009). 



Note. — Cluster positions differ slightly from those reported in Dawson09 due to the use of an updated algorithm for determining cluster centers. 
a photometric redshift 



information for SNe (in the case of clusters A, C, Q, and U). 
We call all visits contingent on the existence of an active SN 
"follow-up" visits (designated by open circles in Fig. 1). 

3. SUPERNOVA SELECTION 

During the survey, our aim was to find as many supernovae 
as possible and find them as early as possible in order to trig- 
ger spectroscopic and NICMOS follow-up. Thus, software 
thresholds for flagging candidates for consideration were set 
very low, and all possible supernovae were carefully consid- 
ered by a human screener. Over the course of the survey, 
thresholds were changed and the roster of people scanning the 
subtractions changed. As a result, the initial candidate selec- 
tion process was inclusive but heterogeneous, and depended 
heavily on human selection. This made it difficult to calculate 
a selection efficiency for the SN candidates selected during 
the survey (listed in Tables 3 and 4 of Dawson09). 

In this section, we select an independent SN candidate sam- 
ple (without regard for the Dawson09 sample) using auto- 
mated selection wherever possible. Although the remainder 
of this paper will focus on cluster SNe, candidates are se- 
lected without regard for cluster membership (which is only 
known from follow-up spectroscopy once the candidate has 
already been found) and we determine SN types for both clus- 
ter and non-cluster SNe. The non-cluster SNe are considered 
further in a second paper deriving the volumetric SN la field 
rate (Barbary et al., in preparation). The automated selec- 
tion consists of initial detection in pairs of subtracted images 
(§3.1; 86 candidates selected), and subsequent requirements 



based on the light curve of each candidate (§3.2; 60 candi- 
dates remaining). The selection efficiency for these two steps 
is later calculated via a Monte Carlo simulation. In §3.3 we 
assign a type (SN la, core-collapse SN, or other) to each of 
the remaining 60 candidates based on all data available (in- 
cluding triggered follow-up observations). For this last step 
we do not calculate an efficiency or completeness. Instead we 
estimate the classification uncertainty of the assigned type for 
each candidate individually. For most candidates the uncer- 
tainty in the type is negligible thanks to ample photometric 
and spectroscopic data. 

3.1. Initial detection 

For the purpose of initially detecting candidates, we use 
only "search" visits (filled circles in Fig. 1) and disregard the 
"follow-up" visits (open circles in Fig. 1). (In the following 
section we will use any available "follow-up" visits to con- 
struct more complete light curves for the candidates discov- 
ered in this section.) We use the MULTlDRlzZLE-combined, 
cosmic ray-rejected, Z850 image from each "search" visit. We 
consider only regions in this image that are covered by three 
or more ZS50 exposures. With less than three exposures, the 
combined images are too heavily contaminated by cosmic 
rays to be practically searchable for SNe. Although there 
are typically four ZS50 exposures, the dither pattern used in 
the survey means that not all regions of the combined image 
have four exposures. The ACS camera is a mosaic of two 
2048 x 4096 pixel CCD chips (1 pixel = 0.05") separated 
by 2.5". The zs^o exposures were dithered to cover this gap, 
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Figure 1. Dates of visits to each cluster. All visits included Z850 exposures 
(usually four). Most visits also included one 2775 exposure. Filled circles in- 
dicate "search" visits (used for finding SNe). Open circles indicate "follow- 
up" visits (contingent on the existence of an active SN candidate). Clusters D, 
N, P, Q and R were re-visited once towards the end of the survey, with addi- 
tional follow-up visits devoted to clusters in which promising SN candidates 
were found (N, Q, R). 



meaning that a5 ;/ wide region in the center of the image and 
2.5" wide regions on either side of the image are only covered 
by two exposures and thus are not searchable. Due to orbital 
constraints, the position angle of HST changes between each 
visit. This means that the unsearchable "gap" region rotates 
over the field between visits, and that the outer parts of the 
field are observed in some visits, but not others (Fig. 2, sec- 
ond row). The regions around bright stars are also considered 
"not searchable" and are similarly masked. 

For each "search" visit to each cluster, we follow these four 
steps: 

1. A reference image is made by combining images from 
other visits to the cluster. All visits that are either 50 or more 
days before the search epoch or 80 or more days after the 
search epoch are included. If there are no epochs outside this 
130 day range, the range is narrowed symmetrically until one 
epoch qualifies. Masked pixels in each visit's image do not 
contribute to the stacked reference image (Fig. 2, third row). 

2. A subtracted image is made by subtracting the stacked 
reference image from the search epoch image. A map of the 
sky noise level in the subtraction is made by considering the 
noise level of the search epoch image and the noise level of 
each reference image contributing to a given region. Any area 
masked in either the search epoch or stacked reference image 
is masked in the subtracted image (Fig. 2, fourth row). 

3. Candidates in the subtraction are identified by soft- 



ware. To be flagged, a candidate must have three contiguous 
pixels with a flux 3.4 times the local sky noise level in the 
subtraction (as determined by the sky noise map above). Once 
flagged, it must fulfill the following four requirements: 

• MULTlDRlzZLE-combined image: A total signal-to- 
noise ratio (including sky and Poisson noise) of 5 or 
more in a 3 pixel radius aperture. 

• MULTlDRlzZLE-combined image: A total signal-to- 
noise ratio of 1.5 or more in a 10 pixel radius aperture. 

• Individual exposures: A signal-to-noise ratio of 1 or 
greater in a 3 pixel radius aperture in three or more in- 
dividual exposures. 

• Individual exposures: A candidate cannot have an indi- 
vidual exposure with a flux more than 20a greater than 
the flux in the lowest flux exposure and a second indi- 
vidual exposure with flux more than 10a greater than 
the flux in the lowest flux exposure. 

The first requirement is designed to eliminate low significance 
detections on bright galaxies. The second requirement helps 
eliminate dipoles on bright galaxy cores caused by slight im- 
age misalignment. The third and fourth requirements are 
aimed at false detections due to cosmic ray coincidence. They 
require the candidate to be detected in most of the exposures 
and allow no more than one exposure to be greatly affected by 
a cosmic ray. On the order of five to ten candidates per sub- 
traction pass all the requirements, resulting in approximately 
1000 candidates automatically flagged across the 155 search 
visits. 

4. Each candidate is evaluated by eye in the subtraction. 

Because the position angle changes between each epoch, the 
orientation of stellar diffraction spikes changes, causing the 
majority of the false detections. These are easy to detect and 
eliminate by eye. Occasionally there are mis- subtractions on 
the cores of bright galaxies that pass the above requirements. 
Only completely unambiguous false detections are eliminated 
in this step. If there is any possibility the candidate is a real 
SN, it is left in the sample for further consideration. 

After carrying out the above four steps for all 155 search 
visit, 86 candidates remain. At this point, candidates have 
been selected based only on information from a single z^^o 
subtraction. Detailed information on each of the 86 candi- 
dates is available from the HST Cluster SN Survey website 28 . 

3.2. Lightcurve Requirements 

The 86 remaining candidates still include a considerable 
number of non-SNe. We wish to trim the sample down as 
much as possible in an automated way, so that we can easily 
calculate the efficiency of our selection. For each candidate, 
we now make three further automated requirements based on 
^775 data (if available) and the shape of the ZS50 light curve. 
The requirements and number of candidates remaining after 
each requirement are summarized in Table 2. 

First, we require that if 2775 data exists for the epoch in 
which the candidate was detected, there be positive flux in a 
2 pixel radius aperture at the candidate location in the 2775 
image. From our SN light curve simulations, we find that vir- 
tually all SNe should pass (near maximum light there is typ- 
ically enough SN flux in the 2775 filter to result in a positive 

28 http://supernova.lbl.gov/2009ClusterSurvey 
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Figure 2. An example of image orientation and searchable regions for cluster ISCS J1432.4+3332. Each column represents an observation of the cluster. The 
first row is the zgso image for that visit. The second row is the part of that image that is searchable. The third row shows the searchable area of the stacked 
reference image used in the subtraction for this visit. The fourth row is the searchable area in the subtraction (the intersection of the second and third rows). 



Table 2 

Light Curve Requirements 



Requirement 


Candidates Remaining 


Before light curve requirements 


86 


Positive ^775 flux (if observed in 2775) 


81 


2cr Detection in surrounding epochs 


73 


If declining, Require two 5<r detections 


60 



total flux, even with large negative sky fluctuations). Mean- 
while, about half of the cosmic rays located far from galaxies 
will fail this test (due to negative sky fluctuations). If there is 
no ^775 data for the detection epoch, this requirement is not 
applied. Even though nearly all SNe are expected to pass, we 
account for any real SNe that would be removed in our Monte 
Carlo simulation. 

Second, we require that the light curve does not rise and 
fall too quickly: if there is a "search" visit less than 60 days 
before the detection visit and also one less than 60 days after 
the detection visit, the candidate must be detected at a 2a level 
in at least one of these two visits. SNe la have light curves 
wide enough to be detected at this level in two epochs spaced 
apart by 60 days. However, cosmic rays in one Z850 image 
are unlikely to be repeated in the same spot in two epochs and 
thus will be removed. This requirement is also included in our 
Monte Carlo simulation. 

The third and final requirement aims to eliminate candi- 
dates that were significantly detected in only the first epoch 
and that then faded from view. Such candidates would not 



have been followed up spectroscopically and it would typi- 
cally be impossible to tell if such candidates were SNe (and if 
so, Type la or core collapse) on the basis of a single detection. 
We chose to eliminate any such candidates and account for 
this elimination in our Monte Carlo simulation, rather than 
dealing with an "untypeable" candidate. Specifically, if a 
candidate is found on the decline (in the first search epoch), 
we require two epochs with 5a detections. For high-redshift 
(z ~ 1) SNe la, this requirement means that the first epoch 
will be at approximately maximum light, and most of the SN 
decline is captured, making it possible to confirm a SN and 
estimate a type. For candidates that are only significantly de- 
tected in the last search epoch, typing is not a problem be- 
cause additional ACS orbits were typically scheduled in order 
to follow such candidates. 

After these requirements 60 candidates remain. The auto- 
matic selection means that we can easily calculate the com- 
pleteness of the selection so far; any real SNe la removed will 
be accounted for in the "effective visibility time" (§4) which 
is calculated using a Monte Carlo simulation. 



3.3. Typing 

We now use all available information about each candi- 
date (spectroscopic confirmation, host galaxy redshift, all 
light curve information, as well as host galaxy luminosity and 
color) to classify each of the 60 remaining candidates as im- 
age artifact, active galactic nucleus (AGN), core-collapse SN 
(SN CC), or SN la. 
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3.3.1. Image artifacts 

Although the automated selections were designed to elim- 
inate image artifacts such as subtraction residuals and cos- 
mic rays, they were made to be somewhat tolerant so that real 
SNe were not eliminated. The result is that some artifacts slip 
through. Candidates located close to the cores of relatively 
bright galaxies that show adjoining negative and positive ar- 
eas in subtractions are likely to be caused by mis-alignment 
between the reference and search image. For such candidates, 
we inspected the full light curve for consistency with the gen- 
eral shape of a SN la light curve. For fourteen of these, the 
light curve is completely inconsistent with that of a SN la. 
Their light curves have either multiple peaks, long flat por- 
tions followed by one or two lower points, and/or 2775 data 
that shows no change. We classify these fourteen candidates 
as subtraction residuals with negligible classification uncer- 
tainty (very unlikely that any are SNe la). 

Candidates where one or two of the four ZS50 exposures 
was clearly affected by a cosmic ray or hot pixel may be false 
detections. These can pass the automated cosmic ray rejec- 
tion when they occur on a galaxy. For two such candidates, 
we used the lack of any change in the 2775 light curve to rule 
out a SN la: fitting SN templates with a range of redshifts 
and extinctions resulted in observed 2775 fluxes too low by 
4cr or more, given the zs^o increase. One other candidate, 
SCP06W50, is less certain. It was discovered in the last visit 
to the cluster, making it difficult to constrain a template light 
curve. There is clearly a hot pixel or cosmic ray in one ZS50 
exposure, but there appears to be some excess flux in the other 
three exposures as well. Also, there is a point-source like de- 
tection in ^775, but offset ~1.2 pixels from the zs^o detection. 
While the 2775 detection may also be a cosmic ray, it is pos- 
sible that this candidate is a SN caught very early. The el- 
liptical "host" galaxy was not observed spectroscopically, but 
we estimate its redshift to be 0.60 < z < 0.85 based on the 
color of ^775 — ^850 = 0.55 and stellar population models of 
Bruzual & Chariot (2003, hereafter BC03). 

Of the 17 total candidates classified as image artifacts, 
SCP06W50 is the only one with significant uncertainty. How- 
ever, this uncertainty does not affect the cluster SN la rate as 
the host galaxy is clearly in the cluster foreground. 

3.3.2. AGN 

Candidates positioned directly on the cores of their host 
galaxies may be AGN. Four such candidates were spec- 
troscopically confirmed as AGN: SCP06L22 (z = 1.369), 
SCP06V6 (z = 0.903) and SCP05X13 (z = 1.642) and 
SCP06U3 (z = 1.534). A fifth candidate, SCP06F3, is spec- 
troscopically consistent with an AGN at z = 1.21, but is less 
certain (see spectroscopy reported in Morokuma et al. 2010). 
SCP06L22, SCP05X13, SCP06U3 and SCP06F3 also have 
light curves that are clearly inconsistent with SNe la (observer 
frame rise times of 100 days or more, or declining phases pre- 
ceding rising phases). Of the "on core" candidates that were 
not observed spectroscopically, five exhibit light curves that 
decline before rising or have rise times of 100 days or more. A 
sixth candidate, SCP06Z51 exhibited slightly varying fluxes 
that could be due to either subtraction residuals or an AGN. 
However, its light curve is clearly inconsistent with a SN la, 
especially considering the apparent size, magnitude and color 
of the host galaxy. Summarizing, there are 1 1 "on-core" can- 
didates certain not to be SNe la. 

Three other "on-core" candidates are also considered 



likely AGN on the basis of their light curves: SCP06Z50, 
SCP06U50 and SCP06D51. These three candidates are 
shown in Figure 3. SCP06Z50 (Fig. 3, top left), has a rise- 
fall behavior in the first three zs^o observations of its light 
curve that could be consistent with a SN la light curve. How- 
ever, given that the host galaxy is likely at z < 1 based on its 
magnitude and color, the SN would be fainter than a normal 
SN la by 1 magnitude or more. Considering the proximity 
to the galaxy core and the additional variability seen in the 
last two observations, SCP06Z50 is most likely an AGN. The 
light curve of candidate SCP06U50 (Fig. 3, top right) also 
exhibits a rise-fall that could be consistent with a supernova 
light curve. However, its host is morphologically elliptical 
and likely at z < 0.7 based on its color. At z < 0.7, a SN la 
would have to be very reddened (E(B — V) > 1) to match 
the color and magnitude of the SCP06U50 light curve. As 
this is very unlikely (considering that the elliptical host likely 
contains little dust), we conclude that SCP06U50 is also most 
likely an AGN. Finally, SCP06D51 (Fig. 3, bottom left) was 
discovered in the last visit, on the core of a spiral galaxy. We 
classify it as an AGN based on the earlier variability in the 
light curve. As these galaxies are all most likely in the clus- 
ter foregrounds, even the small uncertainty in these classifica- 
tions is not a concern for the cluster rate calculation here. 

Note that one of the candidates classified here as a clear 
AGN, SCP06U6, was reported as a SN with unknown red- 
shift by Dawson09, due to the fact that spectroscopy revealed 
no evidence of an AGN. However, it is on the core of a com- 
pact galaxy, and has a clear > 100 day rise in both z^so and 
^775 (Fig. 3, bottom right). While it could possibly be a very 
peculiar SN with a long rise time, what is important for this 
analysis is that it is clearly not a SN la. 

3.3.3. Supernovae 

After removing 17 image artifacts and 14 AGN, 29 candi- 
dates remain (listed in Table 3). One of these is the peculiar 
transient SCP 06F6 (also known as SN SCP06F6) reported by 
Barbary et al. (2009). Various explanations have been con- 
sidered by, e.g., Gansicke et al. (2009), Sokeretal. (2010) 
and Chatzopoulos et al. (2009). It appears that SCP 06F6 
may be a rare type of supernova, with redshift z = 1.189 
(Quimby et al. 201 1; Pastorello et al. 2010). While its precise 
explanation is still uncertain, the important fact for this anal- 
ysis is that SCP 06F6 is clearly not a SN la. 

Note that Table 3 contains 10 fewer candidates than the 
list presented by Dawson09. This is unsurprising; here we 
have intentionally used a stricter selection than in the origi- 
nal search, the source for the Dawson09 sample. Still, after 
finalizing our selection method we checked that there were no 
unexpected discrepancies. Five of the Dawson09 candidates 
(SCP06B4, SCP06U2, SCP06X18, SCP06Q31, SCP06T1) 
fell just below either the detection or signal-to-noise thresh- 
olds in our selection. These were found in the original search 
because detection thresholds were set slightly lower, and be- 
cause the images were sometimes searched in several different 
ways. For example, in the original search SCP06B4 was only 
found by searching an 2775 subtraction. Two Dawson09 can- 
didates (SCP05D55, SCP06Z52) were found too far on the 
decline and failed the light curve requirements (§3.2). Three 
Dawson09 candidates (SCP06X27, SCP06Z13, SCP06Z53) 
were found while searching in "follow-up" visits, which were 
not searched here. SCP06U6 passed all requirements, but is 
classified here as an AGN, as noted above. With the excep- 
tion of SCP06U6, all of these candidates are likely to be su- 
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Figure 3. Images and light curves of four of the 14 candidates classified as AGN. For each candidate, the upper left panel shows the two-color stacked image 
(2775 and Z85o) of the host galaxy, with the position of the transient indicated. The three smaller panels below the stacked image show the reference, new, and 
subtracted images for the discovery visit. The right panel shows the light curve at the SN position (including host galaxy light) in the Z850 (top) and 2775 (bottom) 
bands. The y axes have units of counts per second in a 3 pixel radius aperture. The effective zeropoints are 23.94 and 25.02 for Z850 an d 2775, respectively. The 
discovery visit is indicated with an arrow in the Z85o plot. 



pernovae (mostly core collapse). However, the types of candi- 
dates that did not pass our requirements are not of concern for 
this analysis. Finally, SCP06M50 was not reported in Daw- 
son09, but is classified here as a SN, although a highly uncer- 
tain one (discussed in detail in §3.3.4). 

Thanks to the extensive ground-based spectroscopic follow- 
up campaign, we were able to obtain spectroscopic redshifts 
for 25 of the 29 SNe. The redshift reported in Table 3 is 
derived from the SN host galaxy for all but one candidate 
(SCP06C1) where the redshift is from the SN spectrum it- 
self. Of the 25 candidates with redshifts, eight are in clusters 
and 17 are in the field. Note that this high spectroscopic com- 
pleteness is particularly important for determining the cluster 
or non-cluster status of each SN, which directly affects the 
determination of the cluster SN la rate. The possible clus- 
ter memberships of the four candidates lacking redshifts are 
discussed below. 

We determine the type of each of the 29 supernovae using 
a combination of methods in order to take into account all 
available information for each supernova. This includes (a) 
spectroscopic confirmation, (b) the host galaxy environment, 
and (c) the SN light curve. To qualify the confidence of each 
supernova's type, we rank the type as "secure," "probable," or 
"plausible": 

Secure SN la: Has spectroscopic confirmation or both of the 
following: (1) an early-type host galaxy with no recent 
star formation and (2) a light curve with shape, color 
and magnitude consistent with SNe la and inconsistent 
with other types. 

Probable SN la: Fulfills either the host galaxy requirement 
or the light curve requirement, but not both. 

Plausible SN la: The light curve is indicative of a SN la, but 
there is not enough data to rule out other types. 



Secure SN CC: Has spectroscopic confirmation (note that 
there are no such candidates in this sample). 

Probable SN CC: The light curve is consistent with a core- 
collapse SN and inconsistent with a SN la. 

Plausible SN CC: Has a light curve indicative of a core- 
collapse SN, but not inconsistent with a SN la. 

This ranking system is largely comparable to the "gold," "sil- 
ver," "bronze" ranking system of Strolger et al. (2004), except 
that we do not use their "UV deficit" criterion. This is because 
our data do not include the bluer F606W filter, and because 
SNe la and CC are only distinct in UV flux for a relatively 
small window early in the light curve. Below, we discuss in 
detail the three typing methods used. 

(a) Spectroscopic confirmation: During the survey, seven 
candidates were spectroscopically confirmed as SNe la (Daw- 
son09, Morokumaet al. 2010). These seven (three of which 
are in clusters) are designated with an "a" in the "typing" col- 
umn of Table 3. All seven candidates have a light curve shape, 
absolute magnitude and color consistent with a SN la. Al- 
though the spectroscopic typing by itself has some degree of 
uncertainty, the corroborating evidence from the light curve 
makes these "secure" SNe la. 

(b) Early- type host galaxy: The progenitors of core- 
collapse SNe are massive stars (> 8M ) with main sequence 
lifetimes of < 40 Myr. Thus, core-collapse SNe only occur in 
galaxies with recent star formation. Early-type galaxies, hav- 
ing typically long ceased star formation, overwhelmingly host 
Type la SNe (e.g., Cappellaro et al. 1999; Hamuy et al. 2000). 
In fact, in an extensive literature survey of core-collapse SNe 
reported in early-type hosts, Hakobyan et al. (2008) found 
that only three core-collapse SNe have been recorded in early- 
type hosts, and that the three host galaxies in question had ei- 
ther undergone a recent merger or were actively interacting. 
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Supernovae 










ID 


Nickname 


R.A. (J2000) 


Decl. (J2000) 




SN Type 


Confidence 


Typing 


Cluster Members 
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b 
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la 
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SN SCP06U4 
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la 
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Cluster Membership Uncertain 
















SN SCP06E12 
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02 h 20 m 52 s . 368 


-03° 34' 


13".32 




CC 


plausible 


c 


Not Cluster Members 


















SN SCP06A4 


Aki 


22 h 16 m 01 s . 077 


-17° 37' 


22". 09 


1.193 


la 


probable 


c 


SN SCP06B3 


Isabella 


22 h 05 m 50 s . 402 


-01° 59' 


13". 34 


0.743 


CC 


probable 


c 


SN SCP06C0 


Noa 


12 h 29 m 25 s . 654 
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la 
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CC 


probable 
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SN SCP05D6 


Maggie 
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la 


secure 
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SN SCP06F6 
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a 


SN SCP06F8 
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CC 
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SN SCP06G3 


Brian 


14 h 29 m 28 s . 430 
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la 
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la 
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SN SCP06N33 


Naima 
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1.188 


la 


probable 
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SN SCP05P1 


Gabe 
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la 
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SN SCP05P9 
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0.821 
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1.44 


la 
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SN SCP06Z5 


Adrian 


22 h 35 m 24 s . 966 
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0.623 


la 


secure 


a,c 


Note. — Typing: (a) Spectroscopic confirmation, (b) Host is morphologically early-type, with no signs of recent star formation, (c) Light curve 


shape, color, magnitude consistent 


with type. We do not assig 


n a type for SCP06M50 because there is significant uncertainty that the candidate is a SN at all. 









In all three cases there are independent indicators of recent 
star formation. Therefore, in the cases where the host galaxy 
morphology, photometric color, and spectrum all indicate an 
early-type galaxy with no signs of recent star formation or in- 
teraction, we can be extremely confident that the SN type is la. 
These cases are designated by a "b" in the "typing" column of 
Table 3. We emphasize that in all of these cases, spectroscopy 
reveals no signs of recent star formation and there are no vi- 
sual or morphological signs of interaction. (See Meyers 1 1 for 
detailed studies of these SN host galaxy properties.) 

(c) Light curve: SNe la can be distinguished from most 
common types of SNe CC by some combination of light curve 
shape, color, and absolute magnitude. We compare the light 
curve of each candidate to template light curves for SN la 
and various SN CC subtypes to test if the candidate could 
be a SN la or a SN CC. For candidates lacking both spec- 
troscopic confirmation and an elliptical host galaxy, if there 
is sufficient light curve data to rule out all SN CC subtypes, 
the candidate is considered a "probable" SN la. If SN la can 
be ruled out, it is considered a "probable" SN CC. If neither 
SN la nor SN CC can be ruled out, the candidate is consid- 
ered a "plausible" SN la or SN CC based on how typical 
the candidate's absolute magnitude and/or color would be of 
each type. This approach can be viewed as a qualitative ver- 



sion of the pseudo-Bayesian light curve typing approaches 
of, e.g., Kuznetsova & Connolly (2007); Kuznetsova et al. 
(2008); Poznanski et al. (2007 a,b). SNe classified as "prob- 
able" here would likely have a Bayesian posterior probability 
approaching 1, while "plausible" SNe would have an inter- 
mediate probability (likely between 0.5 and 1.0). We con- 
sciously avoid the full Bayesian typing approach because it 
can obscure large uncertainties in the priors such as lumi- 
nosity distributions, relative rates, light curve shapes, and 
SN subtype fractions. Also, the majority of our candidates 
have more available light curve information than those of 
Kuznetsova et al. (2008) and Poznanski et al. (2007b), mak- 
ing a calculation of precise classification uncertainty less nec- 
essary. In general, classification uncertainty from light curve 
fitting is not a concern for the cluster rate calculation as most 
cluster-member candidates are securely typed using methods 
(a) and/or (b), above. It is more of a concern for the volumet- 
ric field rate calculation based on the non-cluster candidates 
(Barbary et al., in preparation), though the uncertainty in the 
field rate is still dominated by Poisson error. 

For each candidate we fit template light curves for SN la, 
Ibc, II-P, II-L, and Iln. We use absolute magnitude and color 
as a discriminant by limiting the allowed fit ranges accord- 
ing to the known distributions for each subtype. For SN la 
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Table 4 

SN light curve template parameter ranges 



SN type 


Template 


Observed Mb 


E(B 


-v) 


s 


la 


Hsiao 


-17.5- -20.1 


-0.2 


-0.6 


0.6-1.3 


Ibc 


Nugent 


-15.5- -18.5 


-0.1 


-0.5 


1.0 


II-L 


Nugent 


-16.0- -19.0 


-0.1 


-0.5 


1.0 


II-P 


Nugent 


-15.5- -18.0 


-0.1 


-0.5 


1.0 


Iln 


Nugent 


-15.5- -19.1 


-0.1 


-0.5 


1.0 



we start with the spectral time series template of Hsiao et al. 
(2007), while for the core-collapse types we use templates of 
Nugent et al. (2002) 29 . Each spectral time series is redshifted 
to the candidate redshift and warped according to the desired 
color. Observer-frame template light curves are then gener- 
ated by synthetic photometry in the 2775 and z^so filters. The 
magnitude, color, date of maximum light, and galaxy flux in 
^775 and ^850 are allowed to vary to fit the light curve data. 
For the SN la template, the linear timescale or "stretch" (e.g., 
Perlmutter et al. 1997; Guy et al. 2005) is also allowed to vary 
within the range 0.6 < s < 1.3. We constrain the abso- 
lute magnitude for each subtype to the range observed by 
Lietal. (2011); Our allowed range fully encompasses their 
observed luminosity functions (uncorrected for extinction) for 
a magnitude-limited survey for each subtype. We correct from 
their assumed value of Hq = 73 km s - * Mpc -1 to our as- 
sumed value of Ho = 70 km s _1 Mpc -1 and if -correct 
from R to B band. To avoid placing too strict of an up- 
per limit on SN CC brightness, we use the bluest maximum- 
light spectrum available when if-correcting (e.g., for SN Ibc 
we use a bluer spectrum than that of Nugent et al. (2002), as 
bluer SNe Ibc have been observed). The resulting allowed 
Mb range for each subtype is shown in Table 4. Note that 
the range for Ibc does not include ultra-luminous SNe Ic 
(such as those in the luminosity functions of Richardson et al. 
(2002)) as none were discovered by Li et al. (2011). While 
such SNe can mimic a SN la photometrically, the Li et al. 
(2011) results indicate that they are intrinsically rare, and 
even Richardson et al. (2002) show that they make up at most 
~20% of all SNe Ibc. Still, we keep in mind that even can- 
didates compatible only with our SN la template and incom- 
patible with SN CC templates may in fact be ultra-luminous 
SNe Ic, though the probability is low. This is why any candi- 
date typed based on light curve alone has a confidence of at 
most "probable," rather than "secure." The allowed ranges of 
"extinction," E{B — V), are also shown in Table 4. For SN la, 
E(B—V) is the difference in B—V color from the Hsiao et al. 
(2007) template. As the observed distribution of SNe includes 
SNe bluer than this template, SNe la as blue as E(B - V) = 
—0.2 are allowed. Given an E(B — V), the spectral tem- 
plate is warped according to the SALT color law (Guy et al. 
2005), with an effective R B = 2.28 (Kowalski et al. 2008). 
For SN CC templates, extinction as low as E(B — V) = —0.1 
is allowed to reflect the possibility of SNe that are intrinsically 
bluer than the Nugent et al. (2002) templates. Templates are 
then warped using a Cardelli et al. (1989) law with Rb =4.1. 
Extinctions are limited to E(B — V) < 0.5 (implying an ex- 
tinction of Ab ~ 2 magnitudes for SNe CC). 

The light curve template with the largest \ 2 P- value is gen- 
erally taken as the type. We also evaluate each fit by eye to 
check that the best-fit template adequately describes the light 

29 See http://supernova.lbl.gov/~nugent/nugent_templates.html. 



curve. Figure 4 shows the best-fit template for each candidate. 
For candidates typed on the basis of spectroscopic confirma- 
tion or an elliptical host galaxy only the SN la template is 
shown. For candidates typed on the basis of the light curve 
alone, we show both the best-fit SN la and best-fit SN CC 
templates for comparison. The confidence in the best-fit tem- 
plate is either "probable" or "plausible" depending on how 
well other templates fit: If the next-best fit has a P- value that 
is smaller than 10 -3 x Pbest, the best-fit template is consid- 
ered the only acceptable fit and the confidence is "probable." 
If the next-best fit has a P- value larger than 10 -3 x Pbest the 
confidence is "plausible." Finally, note that the photometry 
used here is simple aperture photometry with fixed aperture 
corrections. For SN la cosmology we use color-dependent 
aperture corrections, as described in Suzuki et al. (in prepara- 
tion). 

3.3.4. Comments on individual SN light curves 

Here we comment in greater detail on a selection of individ- 
ual candidates, particularly those with the greatest uncertainty 
in typing. For each candidate, see the corresponding panel of 
Figure 4 for an illustration of the candidate host galaxy and 
light curve. 

SN SCP06E12. We were unable to obtain a host galaxy 
redshift due to the faintness of the host. The color of the 
host galaxy is consistent with the cluster red sequence. The 
candidate light curve is consistent with a SN la at the cluster 
redshift of z = 1.03, but is also consistent with SN II-L at 
z = 1.03. Different SN types provide an acceptable fit over a 
fairly wide range of redshifts. As the SN la template provides 
a good fit with typical parameters, we classify the candidate 
as a "plausible" SN la. However, there is considerable uncer- 
tainty due to the uncertain redshift. 

SN SCP06N32 also lacks a host galaxy redshift. If the 
cluster redshift of z = 1.03 is assumed, the candidate light 
curve is best fit by a SN Ibc template. A SN la template also 
yields an acceptable fit, but requires an unusually red color 
of E(B - V) - 0.6. Given the best-fit s and M B values, 
the candidate would have an unusually large Hubble diagram 
residual of approximately —0.8 magnitudes. If the redshift is 
allowed to float, a SN la template with more typical param- 
eters provides an acceptable fit at z = 1.3. A SN Ibc tem- 
plate still provides a better fit, with the best fit redshift being 
z ~ 0.9. As SN Ibc provides a better fit in both cases, we 
classify this as a "plausible" SN CC. However, there is con- 
siderable uncertainty in both the type and cluster membership 
of this candidate. 

SN SCP06A4. We note that this candidate was observed 
spectroscopically, as reported in Dawson09. While the spec- 
trum was consistent with a SN la, there was not enough evi- 
dence to conclusively assign a type. The host galaxy is mor- 
phologically and photometrically consistent with an early- 
type galaxy, but there is detected [Oil], a possible indication 
of star formation. We therefore rely on light curve typing 
for this candidate, assigning a confidence of "probable" rather 
than "secure." 

SN SCP06G3 has only sparse light curve coverage. The 
best fit template is a SN la with s = 1.3, E(B - V) = 0.3 
and Mb = —18.5, although these parameters are poorly con- 
strained. A large stretch and red color would not be surprising 
given the spiral nature of the host galaxy. It is also consistent 
with a II-L template, although the best fit color is unusually 
blue: E(B -V) = -0.1. Given that SN la yields more "typ- 
ical" fit parameters and that, at z ~ 1 a detected SN is more 
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Figure 4. Images and light curves of the 29 candidates classified as supernovae. For each candidate, the upper left panel shows the 2-color stacked image 
(^775 and Z85o) of the supernova host galaxy, with the SN position indicated. The three smaller panels below the stacked image show the reference, new, and 
subtracted images for the discovery visit. The right panel shows the light curve at the SN position (including host galaxy light) in the zsso (top) and 2775 (bottom) 
bands. The y axes have units of counts per second in a 3 pixel radius aperture. The effective zeropoints are 23.94 and 25.02 for Z850 and 2775, respectively. 
The discovery visit is indicated with an arrow in the Z850 plot. The best-fit SN la template is shown in blue. For cases where the type is SN la based on 
spectroscopic confirmation or host galaxy environment, only the best-fit SN la template is shown, to demonstrate the consistency of the light curve with the 
designation. For cases where the type is based only on the light curve fit, the best-fit core collapse SN template is shown in red. Note that the photometry used 
here is simple aperture photometry with fixed aperture corrections. For SN la cosmology we use color-dependent aperture corrections, as described in Suzuki et 
al. (in preparation). 
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Figure 4. Continued 
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likely to be Type la than II, we classify this as a "plausible" 
Type la, with considerable uncertainty in the type. 

SN SCP06L21 lacks a spectroscopic redshift, but has a dis- 
tinct slowly-declining light curve that rules out a z > 0.6 
SN la light curve. Even the best-fit la template at z = 0.55, 
shown in Fig. 4), is unusually dim (Mb ~ —17.5), making 
it unlikely that the candidate is a lower-redshift SN la. The 
light curve is better fit by a SN II-P template (with the best-fit 
redshift being z = 0.65). We therefore classify the candidate 
as a "probable" SN CC. 

SN SCP06M50 is the most questionable "SN" candidate, 
having no obvious 2775 counterpart to the increase seen in 
^850 • It ma Y m f act be an image artifact or AGN. However, 
it appears to be off the core of the galaxy by ~2 pixels (mak- 
ing AGN a less likely explanation), and shows an increase 
in ^850 flux in two consecutive visits, with no obvious cos- 
mic rays or hot pixels (making an image artifact less likely as 
well). The galaxy is likely to be a cluster member: its color 
and magnitude put it on the cluster red sequence, it is morpho- 
logically early-type, and it is only 19" from the cluster center. 
Under the assumption that the candidate is a supernova and at 
the cluster redshift of z = 0.92, no template provides a good 
fit due to the lack of an 2775 detection and the constraints on 
E(B — V). In particular, a SN la template would require 
E(B — V) > 0.6. (The best-fit template shown in Fig. 4 is 
with E(B - V) = 0.6.) If the redshift is allowed to float, it is 
possible to obtain a good fit at higher redshift (z ~ 1.3), but 
still with E(B — V) > 0.4, regardless of the template type. 
Given the color and early-type morphology of the host galaxy, 
it is unlikely to contain much dust. There is thus no consistent 
picture of this candidate as a SN, and we do not assign a type. 
However, note that the candidate is unlikely to be a cluster 
SN la. 

SN SCP05N10 is the lowest-redshift SN candidate in our 
sample atz = 0.203. Its light curve shape is inconsistent with 
a SN la occurring well before the first observation, and its lu- 
minosity is too low for a SN la with maximum only slightly 
before the first observation. Therefore, we call this a "proba- 
ble" SN CC. For all SN types, the best fit requires maximum 
light to occur well before the first observation, making all fits 
poorly constrained. 

SN SCP06X26 has a tentative redshift of z = 1.44, de- 
rived from a possible [Oil] emission line in its host galaxy. 
Given this redshift, a la template provides an acceptable fit, 
consistent with a typical SN la luminosity and color. How- 
ever, we consider this a "plausible," rather than "probable, " 
SN la, given the uncertain redshift and low signal-to-noise of 
the light curve data. 

3.4. Summary 

In the previous section we addressed the type of all 29 can- 
didates thought to be SNe. However only the cluster-member 
SNe la are of interest for the remainder of this paper. There 
are six "secure" cluster-member SNe la, and two "probable" 
SNe la, for a total of eight. In addition, SCP06E12 is a "plau- 
sible" SN la and may be a cluster member. Two other can- 
didates, SCP06N32 and SCP06M50, cannot be definitively 
ruled out as cluster-member SNe la, but are quite unlikely for 
reasons outlined above. We take eight cluster SNe la as the 
most likely total. It is unlikely that both of the "probable" 
SNe la are in fact SNe CC. We therefore assign a classifica- 
tion error of for each of these, resulting in a lower limit 
of seven cluster-member SNe la. There is a good chance that 
SCP06E12 is a cluster-member SN la, while there is only a 



small chance that SCP06N32 and SCP06M50 are either clus- 
ter SNe la. For these three candidates together, we assign a 
classification error of for an upper limit of nine. Thus, 
8 =b 1 is the total number of observed cluster SNe la. 

4. EFFECTIVE VISIBILITY TIME 

With a systematically selected SN la sample now in hand, 
the cluster SN la rate is given by 



Ej T j L j ' 

where Nsn ia is the total number of SNe la observed in clus- 
ters in the survey, and the denominator is the total effective 
time-luminosity for which the survey is sensitive to SNe Ia in 
clusters. Lj is the luminosity of cluster j visible to the sur- 
vey in a given band. Tj is the "effective visibility time" (also 
known as the "control time") for cluster j. This is the effec- 
tive time for which the survey is sensitive to detecting a SN Ia, 
calculated by integrating the probability of detecting a SN Ia 
as a function of time over the span of the survey. It depends 
on the redshift of the SN Ia to be detected and the dates and 
depths of the survey observations. As each cluster has a dif- 
ferent redshift and different observations, the control time is 
determined separately for each cluster. To calculate a rate per 
stellar mass, Lj is replaced by Mj. 

Equation (1) is for the case where the entire observed area 
for each cluster is observed uniformly, yielding a control time 
T that applies to the entire area. In practice, different areas 
of each cluster may have different observation dates and/or 
depths, resulting in a control time that varies with position. 
This is particularly true for this survey, due to the rotation of 
the observed field between visits and the gap between ACS 
chips. Therefore, we calculate the control time as a function 
of position in each observed field, Tj(x,y). As the cluster 
luminosity is also a function of position, we weight the control 
time at each position by the luminosity at that position. In 
other words, we make the substitution 

t 3 l j^ / T j( x ,y) L j( x ,y)' (2) 

J x,y 

The effective visibility time T at a position (x, y) on the sky 
is given by 

pt — OC 

T(x,y)= / r]*(x,y,t)e(x,y,t)dt. (3) 

J t— — oo 

The integrand here is simply the probability for the survey 
and our selection method to detect (and keep) a SN Ia at the 
cluster redshift that explodes at time t, and position (x,y). 
This probability is split into the probability 77* of detecting 
the supernova and the probability e that the supernova passes 
all "light curve" cuts. As each SN has multiple chances for 
detection, the total probability of detection 77* is a combina- 
tion of the probabilities of detection in each observation. For 
example, if we have two search visits at position (x, y), rj*(t) 
is given by 

V *(t) = r h (t) + (l-r h (t)) V2 (t), (4) 

where r]i (t) is the probability of detecting a SN Ia exploding at 
time t in visit i. In other words, the total probability of finding 
the SN Ia exploding at time t is the probability of finding it 
in visit 1 plus the probability that it was not found in visit 
1 times the probability of finding it in visit 2. This can be 
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generalized to many search visits: The contribution of each 
additional visit to the total probability is the probability of not 
finding the SN in any previous visit times the probability of 
finding the SN in that visit. 

In practice, we calculate T(x,y) in two steps: First, we de- 
termine the probability r] of detecting a new point source in 
a single image as a function of the point source magnitude. 
This is discussed in §4.1. Second, for each (x, y) position in 
the observed area we simulate a variety of SN la light curves 
at the cluster redshift occurring at various times during the 
survey. By considering the dates of the observations made 
during the survey at that specific position, we calculate the 
brightness and significance each simulated SN la would have 
in each ZS50 and 2775 image. We then use our calculation of 
r] as a function of magnitude to convert the observed bright- 
ness into a probability of detecting the simulated SN in each 
observation. The light curve simulation is discussed in §4.2. 
The calculation of cluster luminosities, Lj(x, y), is discussed 
in §5. 

4.1. Detection Efficiency Versus Magnitude 

Here we calculate the probability of detecting a new point 
source as a function of magnitude in a single subtraction. We 
use a Monte Carlo simulation in which artificial point sources 
of various magnitudes are added to each of the individual ex- 
posure images from the survey, before they are combined us- 
ing MultiDrizzle. Starting from the individual exposures 
allows us to test both the efficiency of the MultiDrizzle 
process and our cosmic ray rejection (which uses the flux ob- 
served in the individual exposures). The point sources are 
placed on galaxies in positions that follow the distribution of 
light in each galaxy. Poisson noise is added to each pixel in 
the point source. The altered images are then run through 
the full image reduction and SN detection pipeline used in the 
search, and flagged candidates are compared to the input point 
sources. 

We parameterize the detection efficiency by the ratio of 
point source flux to sky noise. This is a good choice because, 
in most cases, the detection efficiency will depend only on the 
contrast between the point source and the sky noise. However, 
there is an additional dependence on the surface brightness at 
the location of the point source: point sources near the core 
of galaxies will have a lower detection efficiency due to ad- 
ditional Poisson noise from the galaxy. For 0.6 < z < 1.5 
galaxies, we estimate that only ~10% of SNe will fall on re- 
gions where galaxy Poisson noise is greater than the sky noise 
(assuming SNe follow the galaxy light distribution). Still, we 
take this effect into account by splitting our sample of artifi- 
cial point sources into four bins in underlying surface bright- 
ness. The detection efficiency is calculated separately in each 
bin (Fig. 5, top left panel). The first two bins, /i > 22.0 and 
22.0 > fi > 20.6 mag arcsec -2 , correspond to lower surface 
brightnesses where sky noise is dominant. As expected, their 
efficiency curves are very similar. In the third and fourth bins, 
corresponding to higher surface brightness, the Poisson noise 
from the galaxy dominates the sky noise, and the efficiency 
drops as a result. 

For reference, the distribution of sky noise in the sub- 
tractions is shown in Figure 5 (right panel). Nearly all 
the searched area has a sky noise level between 0.006 and 
0.012 counts sec -1 pixel -1 . For a typical value of 0.008, we 
show the corresponding point source z^q magnitude on the 
top axis of the left panel. 

We find that the efficiency curve in each bin is well- 



described by the function 

nM _/?(! + «e- bx )[erf((* - c)M) + 1], x < c 
V( ' ~ \ 1(1 + ae- bx )[evf((x - c)/d 2 ) + 1], x > c ' 

(5) 

where x is the ratio of point source flux to sky noise, and 
a, b, c, d\ and cfo are free parameters. An error function is 
the curve one would expect with a constant cut and Gaussian 
noise, but we find that two different scales (d\ and cfo) in the 
error function, as well as an additional exponential term, are 
necessary to describe the slow rise to r] = 1 at large x. This 
slow rise is due to rarer occurrences, such as cosmic rays co- 
inciding with new point sources. The fitted functions for the 
four bins are plotted in the top left of Figure 5 and reproduced 
in the bottom left of the figure for comparison. We use these 
fitted functions to calculate the effective visibility time in the 
following section. 

4.2. Simulated Lightcurves 

We simulate SN la light curves with a distribution of 
shapes, colors and absolute magnitudes. We use the (origi- 
nal) SALT (Guy et al. 2005) prescription in which the diver- 
sity of SN la light curves is characterized as a two-parameter 
family with an additional intrinsic dispersion in luminosity. 
The two parameters are the linear timescale of the light curve 
("stretch", s) and the B — V color excess, c. For each simu- 
lated SN, s and c are randomly drawn from the distributions 
shown in Figure 6 (solid lines). The stretch distribution is 
based on the observed distribution in passive hosts (Fig. 6, 
left panel, grey histogram) in the first-year Supernova Legacy 
Survey (SNLS) sample (Sullivan et al. 2006). Similarly, the 
color distribution is based on the observed color distribution 
(Fig. 6, right panel, grey histogram) in the first-year SNLS 
sample (Astier et al. 2006). The absolute magnitude of each 
simulated SN is set to 

M B = -19.31 - a(s - 1) + /3c + I (6) 

where —19.31 is the magnitude of an s = 1, c = SN la 
in our assumed cosmology (Astier et al. 2006), a = 1.24, 
P = 2.28 (Kowalski et al. 2008), and / is an added "intrin- 
sic dispersion", randomly drawn from a Gaussian distribution 
centered at zero with a = 0.15 mag. 

We have chosen distributions that represent as accurately 
as possible the full distribution of SNe la occurring in real- 
ity. However, note that the control time is not actually very 
sensitive to the assumed distributions. This is because, for 
the majority of cluster redshifts in the survey, the detection 
efficiency is close to 100% during the time of the survey. Su- 
pernovae would thus have to be significantly less luminous in 
order to change the detection efficiency significantly. In the 
following section §4.3 we quantify the effect on the control 
time arising from varying the assumed SN la properties and 
show that they are sub-dominant compared to the Poisson er- 
ror in the number of SNe observed. All sources of systematic 
errors are also summarized in §6.2. 

To generate the simulated light curves in the observed 
bands, we use the Hsiao et al. (2007) SN la spectral time se- 
ries template. For each simulated SN, the spectral time series 
is warped to match the selected color c and redshifted to the 
cluster restframe. Light curves are generated in the observed 
^775 and ^850 filters using synthetic photometry, and the time 
axis is scaled according to the chosen value of s. 

For each cluster, we calculate T(x,y) in bins of 50 x 
50 pixels (2". 5 x 2". 5). In each bin, we simulate 100 SN 
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Figure 5. Point source detection efficiency in a single subtraction, as a function of the ratio of total point source flux to subtraction noise a 
(counts sec -1 pixel -1 ). The artificial point sources are split into four bins depending on the underlying galaxy surface brightness fi (mag arcsec -2 ) at the 
point source position. The efficiency curve is calculated separately for each bin. In the upper left panel, the four bins are shown, offset for clarity. In the lower 
left panel, the fitted curves are reproduced without offset for comparison. Approximately 72,000 artificial point sources were used in total. The right panel shows 
the distribution of the noise level in the subtractions. The noise level differs by a factor of about two from the deepest to shallowest subtractions searched. 




Figure 6. Left panel: stretch distribution used for simulated SNe (solid line) 
and the stretch distribution of first-year SNLS z < 0.75 SNe in passive 
hosts (Sullivan et al. 2006) (grey histogram). Note that the distribution is not 
changed significantly by cutting the sample at z < 0.6. Therefore we do not 
expect the sample to be significantly Malmquist biased. Right panel: color 
distribution of the first-year SNLS z < 0.6 SNe (Astier et al. 2006) (grey his- 
togram) and the color distribution used for simulated SNe (solid line). The 
dotted lines show alternative color distributions used to assess the possible 
systematic error due to varying amounts of SNe being affected by dust. 



light curves at random positions within the bin. For each sim- 
ulated SN light curve, we shift the light curve in time across 
the entire range of observations, starting with maximum light 
occurring 50 days before the first observation and ending with 
maximum light occurring 50 days after the last observation. 
For each step in time we get the z^^o and 2775 magnitude of 
the SN at every date of observation. From the sky noise maps, 
we know the noise at the position of the simulated SN in ev- 
ery image. Using the curves in Figure 5, we convert the SN 
flux-to-noise ratio to the probability of the SN being detected 
in each ZS50 exposure. (Each simulated SN is also assigned a 
host galaxy surface brightness chosen from a distribution, in 
addition to the randomly selected s, c and I parameters; we 
use the Fig. 5 curve that corresponds to this surface bright- 
ness.) At the same time, we calculate the probability that the 



SN passes our light curve cuts (using both z$5o and 2775 sim- 
ulated magnitudes). Multiplying these two probabilities gives 
the total probability of the simulated SN being included in the 
sample if it peaks at the given date. Integrating the probability 
over time (the entire range of dates) gives the control time for 
each simulated SN. We take the average control time of the 
100 SNe as the value for the given bin. The resulting control 
time map, T(x,y), therefore has a resolution of 2 // .5 x 2 // .5. 
T(x,y) is shown for two example clusters in Figure 7. 



4.3. Effect of Varying SN Properties 

If the real distributions of SN la properties differs signif- 
icantly from those assumed in our simulation, the T(x,y) 
maps we have derived could misrepresent the true efficiency 
of the survey. Above we argued that the effect is likely to 
be small because the detection efficiency is close to 100% for 
most of the survey. Here we quantify the size of the possible 
effect on the control time by varying the assumed distribu- 
tions. 

To first order, changing the assumed distributions of s or 
c or changing the assumed spectral time series will affect 
the detection efficiency by increasing or decreasing the lumi- 
nosity of the simulated SN. To jointly capture these effects, 
we shift the absolute magnitude of the simulated SNe la by 
1q 2 mag and recalculate the control times. To first order, this 
is equivalent to shifting the s distribution by As = 0.2/ a ~ 
0.16 or shifting the c distribution by Ac = 0.2//? ~ 0.09. 
A —0.2 mag shift in absolute magnitude increases the control 
time, decreasing the inferred SN la rate by 6%. A +0.2 mag 
shift decreases the control time, increasing the SN la rate by 
8%. These effects are sub-dominant compared to the Poisson 
error of > 30% in the number of SNe observed. (Sources of 
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Figure 7. Example maps of effective visibility time for clusters ISCS J1432.4+3332 (F) and ISCS J1438. 1+3414 (K). The dot denotes the cluster center and the 
inner and outer circles represent 0.5 Mpc and 1.0 Mpc radius, respectively. The "noise" in these maps is due to the finite number (100) of SNe simulated at each 
position. At lower redshift nearly all simulated SNe are recovered at each position, whereas at higher redshift a sizable fraction of simulated SNe are missed, 
resulting in a higher "noise" level. 



error are summarized in §6.2 and Table 8.) 

For the color distribution, in addition to a simple shift, we 
also quantify the effect of including a smaller or larger frac- 
tion of SNe significantly reddened by dust. In fact, we have 
good reasons to believe that most cluster SNe la will be in 
dust- free environments. A large fraction of the stellar mass 
in the clusters (~ 80%) is contained in red- sequence galax- 
ies expected to have little or no dust. Our spectroscopic and 
photometric analysis (Meyers 1 1) of the red-sequence galaxies 
confirms this expectation. Therefore, for our default c distri- 
bution (Fig. 6, right panel, solid line), we assumed that 20% 
of SNe (those occurring in galaxies not on the red sequence) 
could be affected by dust, and that the extinction of these SNe 
would be distributed according to P{Ay) oc exp(— Ay /0.33) 
[the inferred underlying Ay distribution of the SDSS-II sam- 
ple (Kessler et al. 2009)]. All SNe are assumed to have an in- 
trinsic dispersion in color to match the observed SNLS distri- 
bution at c < 0.3. It might be the case that even fewer SNe are 
affected by dust, or (unlikely) more SNe are affected by dust. 
As extreme examples, we tested two alternative distributions 
(dotted lines in Fig. 6). In the first, we assumed that the SNLS 
sample was complete and characterized the full c distribution, 
with a negligible number of c > 0.4 SNe. This increases 
the control time by only 2%. In the second, we increase the 
fraction of dust-affected SNe from 20% to 50%. Even though 
this alternative distribution includes an additional ~30% more 
reddened SNe (unlikely to be true in reality), the average con- 
trol time is only lower by 9% (increasing the rate by 10%). 
We use these values as the systematic error in the assumed 
dust distribution. 

5. CLUSTER LUMINOSITIES AND MASSES 

In this section, we calculate the total luminosity of each 
cluster and use the luminosity to infer a stellar mass. Only a 
small subset of galaxies in each field have known redshifts, 
making it impossible to cleanly separate cluster galaxies from 
field galaxies. Therefore, we use a "background subtraction" 
method to estimate cluster luminosities statistically: we sum 
the luminosity of all detected galaxies in the field and sub- 
tract the average "background luminosity" in a non-cluster 
field. This approach follows that of Sharon et al. (2007). For 



the blank field, we use the GOODS 30 fields (Giavalisco et al. 
2004) as they have similarly deep or deeper observations in 
both ACS ^775 and zs5o- In §5.1 we describe the galaxy detec- 
tion and photometry method. Simply summing the photome- 
try from the detected galaxies would include most of the total 
cluster light. However, for an unbiased estimate of the total 
light, several small corrections are necessary: We account for 
light in the outskirts of each galaxy (§5.2), and light from faint 
galaxies below the detection threshold (§5.4). These correc- 
tions are on the order of 20% and 5% respectively. In §5.3 we 
convert the observed zs^o flux to a rest- frame 5-band flux. In 
§5.5 we sum the light and subtract background light. In §5.6 
we repeat this calculation limiting ourselves to red- sequence 
and red- sequence early-type subsets of galaxies. Finally, in 
§5.7 we estimate cluster stellar masses based on the cluster 
luminosities and stellar mass-to-light ratios. 

5.1. Galaxy Selection and Photometry 

We use the stacked 2775 and zs^o band images of each 
cluster, which have total exposure times in the range 1060 
- 4450 seconds and 5440 - 16,935 seconds, respectively. 
Galaxy catalogs are created using the method described in de- 
tail by Meyersll: We run SEXTRACTOR (Bertin & Arnouts 
1996) in dual-image mode using the zsso image for detection, 
and use a two-pass Cold/Hot method (Rix et al. 2004) to op- 
timally de-blend galaxies. We remove stars from the catalog 
based on the CLASS_STAR and FLUX.RADIUS parameters 
from the z^so image. 

It is notoriously difficult to determine accurate total fluxes 
for extended sources. However, as we are only concerned 
with the summed flux of many galaxies, it is not important 
that the estimate be accurate for each individual galaxy, only 
that the estimate is unbiased in the aggregate. We use the 
SEXTRACTOR MAG_AUTO photometry (which gives the to- 
tal flux within a flexible elliptical aperture) and apply a correc- 
tion determined using the Monte Carlo simulation described 
below. In order to make the aperture correction as small as 
possible, we use a relatively large "Kron factor" of 5.0, mean- 

30 Based on observations made with the NASA/ESA Hubble Space Tele- 
scope. The observations are associated with programs GO-9425, GO-9583 
and GO-10189 
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Figure 8. Percentage of simulated galaxies recovered by S Extractor as 
a function of total galaxy zsso magnitude for simulated galaxies placed on 
cluster fields (black circles) and GOODS fields (grey squares). The detection 
efficiency drops to 80% at zsso = 24.72 for cluster fields (vertical line). We 
discard all galaxies dimmer than this value. 

ing that the MAG_AUTO aperture is scaled to 5.0 times the 
Kron radius of the galaxy. MAG_AUTO is only used to de- 
termine ^850 magnitudes; 2775 — ZS50 colors are determined 
using PSF matching and a smaller aperture, as described in 
Meyersll. 

5.2. Galaxy Detection Completeness and Magnitude Bias 

To count all the flux in all cluster galaxies, we must make 
two corrections: (1) add the galaxy light outside of the 
MAG_AUTO aperture, and (2) add the luminosity of all clus- 
ter galaxies below the detection threshold of our galaxy cat- 
alog. We use a Monte Carlo simulation of galaxies placed 
on our real survey data to determine both the detection effi- 
ciency as a function of galaxy magnitude, and the fraction of 
galaxy light inside the MAG_AUTO aperture. Each simulated 
galaxy has a Sersic (1968) profile, with the Sersic index n 
simply selected from a flat distribution ranging from n = 0.7 
to n = 4.5, and the minor to major axis ratio q selected from 
a flat distribution ranging from q = 0.3 to q = 1. The dis- 
tribution of galaxy angular sizes will also affect the results. 
For guidance on the size of the galaxies of concern (namely, 
those at z > 0.9) we turned to the subsample of the 672 galax- 
ies having spectroscopic redshifts 0.85 < z < 1.6. These 672 
galaxies were all fit with GALFIT (Peng et al. 2002), which fits 
a value for r e . Based on the distribution of r e as a function of 
magnitude for these galaxies, we chose r e for each simulated 
galaxy (based on its magnitude). A total of 15000 and 12000 
simulated galaxies were placed on cluster and GOODS fields 
respectively. 

The detection efficiency as a function of galaxy magnitude 
is shown in Figure 8. For the average of all cluster fields, 
the detection efficiency drops to 80% at z^so = 24.72. We 
use this magnitude as a cutoff in our selection, discarding 
all galaxies dimmer than this magnitude. We later correct 
total cluster luminosities for the uncounted light from these 
galaxies by using an assumed cluster luminosity function. In 
reality, the detection efficiency varies slightly from field to 
field (and even within a field) due to exposure time variations. 
However, to first order, the variation is accounted for by us- 
ing the average efficiency in all fields. In addition, the total 
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Figure 9. Galaxy MAG.AUTO aperture correction as a function of galaxy 
magnitude. Black circles: Average correction for the full distribution of 
galaxies simulated, including all Sersic indices n. The black line is a fit to 
these points and is the relation we use. Note that it is not extrapolated beyond 
the range shown. To illustrate the effect of n on the aperture correction, we 
plot the aperture correction for subsets of galaxies with different Sersic in- 
dices (Grey squares and triangles). Galaxies with larger Sersic indices have 
a larger aperture correction. 

luminosity of z^so > 24.72 cluster galaxies is expected to be 
small (as we show below), so slight changes in the cutoff will 
have a negligible effect on the total luminosity. 

For each simulated galaxy, we determine the difference 
(AM) between the MAG_AUTO magnitude and the true 
total magnitude. Binning the simulated galaxies by their 
MAG_AUTO magnitude, we derive a relation between AM 
and the galaxy brightness (Fig. 9, black circles). AM gen- 
erally increases with galaxy magnitude because the outskirts 
of dimmer galaxies are increasingly buried in noise, caus- 
ing S Extractor to underestimate the true extent of the 
galaxy, and thereby underestimate the Kron radius, resulting 
in a smaller MAG_AUTO aperture. We find that the relation 
is well-fit by a second-order polynomial (Fig. 9, thick black 
line), given by 



AM = 0.238 + 0.081 (M M ag.auto ~ 23) 
+0.009(M M AG_At/TO-23) 2 . 



(7) 



We use this to correct the magnitude of each detected galaxy. 
Note that the correction is not extrapolated beyond the fitted 
range shown. 

Because we cannot reliably determine r e or the Sersic in- 
dex n for each galaxy, we rely on the simulated distribution 
of r e and n to accurately represent the true distributions. (The 
black circles in Fig. 9 include all simulated galaxies.) We 
have based our distribution of r e on actual galaxies, but n is 
less well-known. To estimate the effect of varying the n dis- 
tribution, we show AM for subsets of the simulated galaxies, 
divided by Sersic index (Fig. 9, grey points and lines). AM 
increases with Sersic index, because a larger Sersic index im- 
plies a larger fraction of light in the outskirts of the galaxy, 
under the detection threshold. This leads to a smaller esti- 
mate of the Kron radius, and a smaller MAG_AUTO aperture. 
If, instead of the flat 1 < n < 4 distribution used, all galax- 
ies had 1 < n < 2, the aperture correction would be lower 
by approximately 0.10 magnitudes. If instead all galaxies had 
3 < n < 4, the correction would be higher by approximately 



The SN la Rate in High-Redshift Galaxy Clusters 



19 



0.07 magnitudes. We use 0.07 mag as the systematic uncer- 
tainty in the aperture correction. (All systematic uncertainties 
are summarized in §6.2 and Table 8.) 

5.3. K -Corrections 

We use a if -correction based on the BC03 stellar popula- 
tion spectral models to convert the observed zgso magnitude 
to a rest-frame B magnitude for each cluster. Rather than 
using a single if -correction for all the light in each cluster, 
we apply a if -correction to each galaxy magnitude based on 
its ^775 — ^850 color. For each cluster's redshift, we deter- 
mine the relation between if -correction (Mb (rest) — £85o) 
and ^775 — ^850 color, using BC03 spectra with initial metal- 
licities in the range 0.004 < Z < 0.05 and ages in the range 
1 x 10 8 — 5 x 10 9 yr. For most cluster redshifts in our sam- 
ple, all of the spectra over this wide range fall along the same 
line in if -correction versus color, meaning that the color de- 
termines the if -correction, regardless of the metallicity or age 
assumed. The dispersion of the models about the best-fit line 
is < 0.03 mag at redshifts < 1.1 and > 1.4, and reaches 
its largest value of 0.09 mag at z = 1.26. We calculate the 
if -correction for each galaxy using this best-fit relation, ef- 
fectively assuming that every galaxy is at the cluster redshift. 
This results in an incorrect luminosity for non-cluster mem- 
ber galaxies, but this is accounted for by performing the same 
if -correction on the galaxies in the GOODS fields prior to 
subtracting their luminosity. 

5.4. Luminosity function correction 

We estimate the total luminosity of all galaxies below the 
detection limit of zs^o = 24.72 using a Schechter (1976) lu- 
minosity function, which gives the number of galaxies in the 
luminosity interval [L,L + dL] in a given sample, 



$(L)dL = $*(L/L*) a e- L/L *d(L/L*). 



(8) 



<£* is a normalization, L* is a characteristic galaxy luminos- 
ity, and a is a unit-less constant. The ratio of total to observed 
luminosity is then 



J °°L^L)dL 



(9) 



and we multiply each observed cluster luminosity by C to get 
the total luminosity. 

We assume values for L* and a determined in other studies 
and use our data to perform a rough consistency check. For 
a, studies have shown that the value does not evolve much 
from low redshift, at least for redder galaxies. Analyzing only 
red galaxies in 28 clusters spanning < z < 1.3, Andreon 
(2008) find a = -0.91 ± 0.06 (rest-frame F-band) with 
no discernible trend in redshift (see also Andreon 2006b,c). 
From five intermediate-redshift clusters (0.54 < z < 0.9), 
Crawford et al. (2009) find a somewhat flatter faint-end slope 
a ~ —0.6 (rest- frame 5-band) for the red-sequence lumi- 
nosity function. Looking at the full luminosity function, 
Goto et al. (2005) find a = -0.82 ± 0.10 in one cluster at 
z = 0.83 (rest- frame 5-band), compared to a = — 1.00±0.06 
in 204 low-redshift clusters (rest-frame g-band) (Goto et al. 
2002). In redder bands, Strazzullo et al. (2006) find a ~ -1 
for three clusters at redshifts 1.11 < z < 1.27 (in approxi- 
mately rest- frame z band). Summarizing, most studies find a 
value consistent with a ~ —0.9, and we assume this value in 
computing C. 



Values for M* are also reported in most of the above- 
mentioned studies. Studies of red galaxies find that the vari- 
ation of M* with redshift is consistent with passive evolu- 
tion, with M* decreasing towards higher redshifts (Andreon 
2006c; Crawford et al. 2009). Crawford et al. (2009) find 
M B = —21.1 and M B ~ -21.3 (with errors of approxi- 
mately a half magnitude) for two clusters at redshifts 0.75 and 
0.83. if -correcting from the observed [3.6] -band, Andreon 
(2006c) find M B ~ -21.7 at z ~ 1.1, with approximately 
0.5 magnitudes of evolution between z = 0.3 and z = 1.1. 
At lower redshift (considering all galaxies) Goto et al. (2002) 



find M? 



-21.6, compared to M B ~ —21.0 for one cluster 



at z = 0.83 (Goto et al. 2005). On the basis of these measure- 



ments, we assume a value of Mt 



-21.7. 



We have checked our assumed M B and a for consistency 
with our data. With the set of spectroscopically-confirmed 
cluster galaxies from our clusters at z < 1.2, we confirmed 
that the bright end of the luminosity function is consistent 
with M B = —21.7, and strongly inconsistent with values out- 
side the range M B = —2 1.7 ±0.5. We also determined the lu- 
minosity function using a statistical subtraction of the "back- 
ground" luminosity function from the GOODS fields, finding 
excellent agreement with the assumed M B and a values over 
the range -24 < M B < -19.8 (M B = -19.8 corresponds 
to the detection limit in the highest-redshift clusters). 

For each cluster, we calculate C in the observer frame, con- 
verting M B = —21.7 to the observed zsso band, using the 
cluster redshift and a if -correction based on a passive galaxy 
template. In Table 5 we report the value Zg 50 and the re- 
sulting correction C for each cluster. The correction is less 
than 5% for the majority of clusters, rising to a maximum of 
14% for the highest-redshift cluster. Because the correction 
is so small, varying the assumed values of M B and a does 
not have a large effect on the total luminosity. Varying M B 
by ±0.5 mag (a larger range than that allowed by our data) 
changes the average correction by only ~^_\%. Varying a by 
±0.2 changes the average correction by ~^%. We conserva- 
tively take ^3°% (the full range when varying both concur- 
rently) as the systematic uncertainty in luminosity from the 
faint-end correction (summarized in §6.2). 

5.5. Cluster Luminosities and Aggregate Cluster Profile 

For each cluster we sum the if -corrected 5-band luminos- 
ity of all galaxies brighter than the detection limit zs5o = 
24.72. To reduce noise, we discard galaxies that are clearly 
too bright to be cluster members. In clusters with a central 
dominant (cD) galaxy or dominant (but not central) bright- 
est cluster galaxy(BCG), the bright cutoff magnitude is set to 
the magnitude of the cD galaxy or BCG. In clusters lacking 
a clearly dominant galaxy, we conservatively set the cutoff 
based on the absolute magnitude of the most luminous cD 
galaxy in any cluster, M B = —23.42 (from cluster XMMU 
J2235.3— 2557). The bright cutoff magnitude in the observer 
frame, £g5Q ght , is listed for each cluster in Table 5. Because 
the bright cutoff is chosen so conservatively, we expect that 
no cluster galaxies are discarded. The effect of being overly 
conservative is only to add noise, and this is captured in the 
statistical uncertainty described below. 

For each cluster we apply the same selection criteria and 
if -corrections to the GOODS fields to determine the "back- 
ground" specific to that cluster. The error in the luminos- 
ity comes from the error in this background determination, 
which we estimate in the following way: We select 30 non- 
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Table 5 

Bright cutoff magnitudes and luminosity function parameters 



ID 


z 


Cutoff from 


bright 
850 


7* 

850 


C 


A 


1.46 


Max cD 


21.09 


22.80 


1.143 


B 


1.12 


cD 


20.11 


21.38 


1.033 


C 


0.97 


cD 


19.87 


20.79 


1.018 


D 


1.02 


BCG 


20.13 


20.95 


1.021 


E 


1.03 


cD 


19.40 


20.99 


1.022 


F 


1.11 


Max cD 


19.63 


21.34 


1.031 


G 


1.26 


BCG 


20.34 


22.04 


1.064 


H 


1.24 


BCG 


20.33 


21.95 


1.058 


I 


1.34 


Max cD 


20.66 


22.37 


1.092 


J 


1.37 


Max cD 


20.77 


22.50 


1.104 


K 


1.41 


Max cD 


20.92 


22.65 


1.122 


L 


1.37 


Max cD 


20.77 


22.50 


1.104 


M 


0.90 


Max cD 


18.69 


20.53 


1.014 


N 


1.03 


BCG 


20.22 


20.99 


1.022 


P 


1.1 


Max cD 


19.58 


21.29 


1.030 


Q 


0.95 


cD 


20.01 


20.66 


1.015 


R 


1.22 


Max cD 


20.15 


21.86 


1.054 


S 


1.07 


Max cD 


19.44 


21.16 


1.026 


T 


0.97 


Max cD 


19.00 


20.75 


1.017 


U 


1.04 


Max cD 


19.31 


21.04 


1.022 


V 


0.90 


cD 


18.89 


20.49 


1.013 


w 


1.26 


Max cD 


20.33 


22.04 


1.064 


X 


1.10 


Max cD 


19.58 


21.34 


1.031 


Y 


1.24 


cD 


20.29 


21.90 


1.056 


Z 


1.39 


cD 


20.85 


22.58 


1.112 



Note. — "Cutoff from" refers to how Zg5 sht is determined. "cD": magnitude of 
visually central dominant galaxy. "BCG": magnitude of visually classified brightest 
cluster elliptical (but not central) galaxy. "Max cD": Cluster does not have obvious cD 
galaxy or clear BCG. In this case, z^ 1 ^ is K-corrected from M B = —23.42, the 
absolute magnitude of the brightest cD galaxy in the entire sample. 



connected circular regions (15 in each of GOODS North and 
South) of radius 1.4', similar to the size of the cluster fields. 
We determine the luminosity density in each of these fields. 
The average is taken as the background luminosity for the 
cluster, and the standard deviation (typically 15 - 20 % of the 
average) is taken as the error in this "background" luminosity 
due to variations between fields. 

We have implicitly assumed that the GOODS average ac- 
curately represents the cosmic average. GOODS incorporates 
only two widely separated fields. As a result, the average lu- 
minosity density may differ from the cosmic average due to 
variations in large scale structure. As a rough estimate of the 
cosmic variance, we compare the two GOODS fields. The 
average luminosity density of the GOODS-North regions is 
consistently higher than that of the GOODS-South regions by 
15 - 20%. This means that the "standard deviation" of these 
two samples of large scale structure is ~8%. We checked 
this using the cosmic variance calculator made available by 
Trenti & Stiavelli (2008) 31 . The expected cosmic variance in 
galaxy number counts in the redshift window 0.7 < z < 1.7 
for one GOODS field is approximately ~6%, in good agree- 
ment with our naive estimate. Conservatively, we take 8% as 
the cosmic variance for one GOODS field. For the average of 
the North and South fields, this implies a cosmic variance of 

8%/y/2 ~ 6%. 

One might be additionally concerned that the "background" 
in the cluster fields is biased higher than the cosmic average 
because clusters form in regions of large-scale overdensities. 
However, each cluster field is a "pencil-beam" galaxy survey, 
so the vast majority of non-cluster galaxies will not be as- 
sociated with the high-density region in which each cluster 

3 1 http://casa.colorado.edu/~ trenti/Cosmic Variance.html 



formed. 

Ideally one would measure a two-dimensional luminosity 
density, L(x,y), for each cluster, as in Equation (2). How- 
ever, the large background makes this difficult. For our pur- 
pose (which is to account for variations in control time with 
radius), it is sufficient to assume the clusters have a circu- 
larly symmetric luminosity distribution, L(r). For each clus- 
ter, we sum the total luminosity in annuli of width 0.1 Mpc. 
For nearly all clusters there is a clear overdensity relative to 
the background out to r ~ 0.3 Mpc. Beyond 0.3 Mpc, the lu- 
minosity measurement is dominated by background noise for 
most clusters. This might appear to be a problem; we wish 
to characterize the cluster luminosities out to r > 0.7 Mpc, 
the area over which we searched for SNe. In fact, it is only 
necessary to accurately measure the average luminosity pro- 
file over the full area (the denominator of Eq. 1 is the sum of 
the cluster luminosities, weighted by control time). Averag- 
ing all 25 clusters, there is a significant measurement of the 
luminosity profile out to > 0.5 Mpc (Fig. 10, left panels), and 
the average cluster luminosity within r < 0.6 Mpc has an 
error of 12% (statistical only) and ~ 20% (statistical + cos- 
mic variance), below the Poisson error in the number of SNe 
detected. 

Beyond r < 0.6 Mpc, the control time is generally small 
(that is, there are few observations covering the outskirts of 
the clusters) and the cluster luminosity density is low, mean- 
ing that these regions will not contribute greatly to the rate 
measurement. Still, we include these regions in our rate cal- 
culation, using the entirely reasonable prior that the luminos- 
ity density is decreasing with radius past r < 0.6 Mpc. How 
rapidly the luminosity density decreases will not have a sig- 
nificant impact on the result, but as a convenient analytic de- 
scription we fit a /3-model of the form 

L (r) = r, T~T~ \oVr ( 10 ) 

W (1 + (r/r core ) 2 )P 

over the range r < 0.6 Mpc and apply this function at 
r > 0.6 Mpc. The data are well-fit by this model, with best- 
fit parameters r core = 0.074 Mpc and f3 = 0.91. Varying this 
model luminosity by AE = ±20% (easily enclosing the al- 
lowed range of L(r)) only changes our results by ±4%. This 
and other systematic uncertainties are summarized in Table 8. 

5.6. Galaxy subsets 

In addition to measuring the total luminosity of all galaxies 
in the clusters, we also measure the total luminosity of only 
red- sequence galaxies and the total luminosity of only red- 
sequence, morphologically early-type galaxies. These mea- 
surements enable us to compute the cluster SN la rate specifi- 
cally in these galaxy subsets. For the red- sequence-only mea- 
surement we follow the same procedure as above, but elim- 
inate from the analysis all galaxies with 2775 — zs^o colors 
more than 0.2 mag from their respective cluster red sequences 
(galaxy colors and cluster red sequences are determined as 
in Meyers 1 1). For the red-sequence early-type measurement, 
we make the same requirement in color, and additionally use 
the quantitative morphology requirements of Meyers 1 1 . Mey- 
ers 11 use two parameters, asymmetry and Gini coefficient, to 
automatically divide galaxies into early- and late-type sub- 
sets. Here we require the asymmetry to be < 0.10 and the 
Gini coefficient to be > 0.40. We also require the galaxies 
to be 2:850 < 24 as the asymmetry and Gini coefficient are 
somewhat less reliable at fainter magnitudes. 
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Figure 10. Average luminosity profile of the 25 clusters. Top row: Average luminosity density in the cluster fields in annuli of width 0.1 Mpc extending out from 
the cluster center. The grey line and shaded region show the estimated "background" luminosity in each annulus and the error on that background, respectively. 
The darker grey region is the statistical-only error, while the light grey is the statistical + cosmic variance error, added in quadrature. Bottom row: The total 
enclosed luminosity as a function of radius, derived by subtracting the background from the total luminosity density in each bin in the top row plot. The left plots 
include galaxies of all colors and morphologies, while the center plots include only galaxies with 2775 — z^bQ colors within ±0.2 mag of the red sequence in 
their respective clusters. The right plots include only galaxies that satisfy the color requirement and also have zsso < 24 and are morphologically early type. By 
excluding bluer galaxies (center and right plots) the background (and error) is reduced dramatically. 



Table 6 



Average cluster luminosities within r < 0.6 Mpc 


Cluster subset 


^clusters 


z 


All galaxies (10 12 L O)jB ) 


RS galaxies (10 12 L O)jB ) 


RSE galaxies (1O 12 L 0jjB ) 


X-ray discovered 


9 


1.20 


2.86 ±0.54 ±0.45 


2.42 ±0.16 ±0.05 


1.47 ±0.12 ±0.02 


IR-Spitzer discovered 


7 


1.30 


2.85 ±0.70 ±0.52 


1.83 ±0.24 ±0.07 


0.96 ±0.16 ±0.03 


Optical discovered 


9 


1.00 


1.99 ±0.37 ±0.32 


1.75 ±0.08 ±0.03 


1.29 ±0.06 ±0.01 


z < 1.2 


14 


1.02 


2.14 ±0.31 ±0.33 


1.79 ±0.07 ±0.03 


1.28 ±0.05 ±0.01 


z > 1.2 


11 


1.32 


3.06 ±0.58 ±0.54 


2.31 ±0.19 ±0.07 


1.23 ±0.14 ±0.04 


All Clusters 


25 


1.15 


2.54 ±0.31 ±0.42 


2.02 ±0.09 ±0.05 


1.26 ±0.07 ±0.02 



Note. — "RS": galaxies within ±0.2 mag of the cluster red sequence. "RSE": galaxies fulfilling the "RS" requirement, and also zsso < 24, and morphologically early-type. The 
first and second confidence intervals are the statistical error and cosmic variance error, respectively. These luminosities do not include the faint-galaxy correction C. 



The luminosity profiles for these two subsets are shown in 
the center and right columns of Figure 10. The profiles are 
broadly consistent with the profile of the full cluster lumi- 
nosity (left column), but the "subset" profiles are much better 
measured. This is because by excluding bluer galaxies, we 
have eliminated much of the background while still retaining 
the majority of cluster galaxies. The red-sequence subset con- 
tains 77% of the luminosity of the full cluster within 0.6 Mpc 
(Table 6). The red-sequence early-type subset has 62% of the 
light contained in the red-sequence subset. However, keep in 
mind that in the early-type subset we have excluded ^50 > 24 
galaxies, whereas they are included in the red- sequence sub- 
set: In fact 68% of Z850 > 24 red-sequence galaxies pass the 
"early-type" morphology requirements. 

Note that our definition of "red- sequence" here is a rela- 
tively simple one. It is sufficient to select a subsample of 
"more red" galaxies for the purpose of looking for a depen- 



dence of the SN rate with galaxy color within the cluster. 
However, for measuring the red fraction in clusters (e.g., the 
Butcher-Oemler effect Butcher & Oemler 1978, 1984), defin- 
ing red- sequences with a constant color width for all redshifts 
is not ideal (Andreon 2006a). The luminosity content of the 
subsets are reported above only to give the relative size of 
each sample; a full analysis of the cluster content is beyond 
the scope of this paper. 

5.7. Stellar Mass -to -Light Ratio 

To compare SN rates in clusters of different ages, rate mea- 
surements must be normalized by stellar mass rather than stel- 
lar luminosity because luminosity changes as stars age. To 
convert our luminosity measurements to mass measurements 
we use a mass-to-light (M/L) ratio based on a stellar evo- 
lution model. There are several available models in the lit- 
erature. The choice of stellar tracks, metallicity, star forma- 
tion history, and in particular the assumed IMF, will all af- 
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feet the derived Mj L ratio to some extent. For the purpose 
of measuring the change in rate with redshift, it is important 
to use a consistent model and assumptions for determining 
the M/L ratio for all rate measurements. That is, we are 
most concerned that the model accurately captures the evo- 
lution of stellar luminosity over the redshift range of interest 
(0 < z < 1.46), and less concerned about the overall nor- 
malization of the M/L ratio. To that end, for our main result 
we will use a model and assumptions that match as closely 
as possible those used for the Mj L ratio in low-redshift clus- 
ter rate measurements. As we also give results normalized by 
luminosity, those wishing to use a different M/L ratio can 
easily do so. Finally, note that the initial stellar mass formed 
is the quantity of interest for normalizing rate measurements. 
However, as most rate measurements and M/L ratios have 
been reported in terms of current mass, we give our results 
in these units and simply note the difference between current 
and initial mass for the purpose of comparing rate measure- 
ments. Thus, in the following paragraphs M refers to current 
stellar mass. 

5.7. 1 . M/ L ratio in low-redshift cluster rate measurements 

The lower-redshift cluster rate studies of Sharon et al. 
(2007), Sharon et al. (2010), and by extension, Dilday et al. 
(2010) have used the relations between Mj L ratio and galaxy 
color derived by Bell et al. (2003, hereafter Bell03). For ex- 
ample, Sharon et al. (2007) use the relation \og 10 (M/L z ) = 
-0.052 + 0.923(r - i) and Sharon et al. (2010) use 
\og 10 (M/L 9 ) = — 0.499+1. 519(# — r), where M, L z <mdL 9 
are in solar units. In order to use a consistent model, it is im- 
portant to recognize how these relations were derived. Bell03 
fit a grid of PEGASE2 (Fioc & Rocca-Volmerange 1997) syn- 
thetic galaxy spectral energy distributions (SEDs) to actual 
ugrizK photometry of low-redshift galaxies. The grid cov- 
ers a range of metallicities and star formation histories. The 
star formation histories have exponentially-decreasing or - 
increasing star formation rates, and assume that star forma- 
tion commenced at z = 4. For each galaxy, the M/L ratio is 
that of the best-fit synthetic galaxy SED, consistently evolved 
to z = 0. Bell03 use a "diet" Salpeter (1955) IMF (following 
Bell & de Jong 2001). This IMF is defined as having the same 
colors and luminosity as a Salpeter IMF, but with a total mass 
30% lower. The difference in mass is attributed to a smaller 
number of faint low-mass stars relative to a Salpeter IMF. 
These stars do not contribute significantly to the luminosity 
of the Salpeter IMF. The diet Salpeter IMF results in M/L 
ratios 30% lower at a given color than a normal Salpeter IMF. 
Note that because Bell03 simply take the M/L ratio from the 
best-fit synthetic SED of each galaxy, the Bell03 relations will 
generally fall within the grid of M/L versus color covered by 
the synthetic galaxy SEDs. 

5.7.2. M/L ratio at 0.9 < z < 1.46 

Ideally, for consistency with Sharon et al. (2007), 
Sharon et al. (2010) and Dilday et al. (2010), we 
would simply use the Bell03 relation for u — g color, 
which most closely matches our observed color: 
\og 10 (M/L g ) = -0.221 + 0A85(u - g). However, 
the Bell03 relations are based on ugrizK photometry of 
low-redshift galaxies, corrected for evolution to z = 0. As 
such, they are specific to z = and not directly applicable 
at high redshift. A stellar population passively evolving 
from age a few Gyr (at z <~ 1) to > 10 Gyr (at z = 0) 



will dim significantly while only growing slightly redder 
(see, e.g. BC03), in a manner that does not follow the 
Bell03 relations. To estimate the effect of evolution from 
their z = relation to higher redshift, we make a similar 
grid of PEGASE2-generated SEDs with the same formation 
redshift, metallicities, IMF, and star formation histories. As 
expected, when evaluated at z = 0, the M/L ratios of this 
grid are consistent with the Bell03 relation (Fig. 11, left 
panel, upper grid of black points). Evaluating the SEDs at 
higher redshifts, we find that the M/L ratios are well fit by 
a relation with the same slope, but smaller normalization. 
For example, at z = 1.2, the best-fit offset from the z = 
relation is —0.36 dex (Fig. 11, left panel, dashed line). At the 
extremes of the redshift range of interest, the best fit offset is 
-0.26 dex (z = 0.9) and -0.44 dex (z = 1.46). We therefore 
use a M/L ratio of 

loe (MIL ) - ( "°- 48 + °' 485 ^ -9)' z = °' 9 

(11) 

and linearly interpolate for intermediate redshifts. Another 
way to view Equation (11) is that, independently of the rela- 
tion at z = 0, we have fit a linear relation to the PEGASE2 
SEDs at the redshift of each cluster, assuming a slope consis- 
tent withBelKB. 

Using Equation (11) we calculate mass on a galaxy-by- 
galaxy basis: we if -correct the observed 2775 and ZS50 magni- 
tude to rest- frame SDSS u and g magnitudes using the method 
discussed in §5.3, and obtain the M/L ratio from the u — g 
color. In all, 66% of the clusters' luminosity is from galax- 
ies with color in the range 1.3 < u — g < 1.7, 27% of the 
luminosity is distributed roughly equally between galaxies in 
the range 0.6 < u — g < 1.3, and the remainder is in red- 
der galaxies with u — g > 1.7. Thus, while there is a clear 
presence of bluer cluster galaxies, the majority of the clusters 
luminosity is confined to a narrow range in color. This nar- 
row color range means that changes in the assumed slope of 
Equation (11) will not have a large effect on the resulting total 
mass. 

The cumulative M/L ratio (the ratio of the total mass of 
all 25 clusters to the total luminosity of all 25 clusters) is 
M/Lg = 1.25 (see Table 7, "denom"). For red-sequence 
galaxies only, the ratio is higher (M/L g = 1.38) due to the 
exclusion of bluer galaxies with a lower inferred M/L ratio. 

5.7.3. M/L ratio uncertainty 

As noted above, we are primarily concerned with the accu- 
racy of the evolution in the stellar mass and luminosity over 
the range < z < 1.46, rather than the accuracy of the abso- 
lute M/L ratio. As a cross-check of the M/L ratio evolution, 
we have compared the above results (using PEGASE2) to the 
results obtained with the BC03 SEDs. We use the standard 
Padova 1994 evolution and the same star formation histories 
as above. In terms of evolution offset from z = to z ~ 1.2, 
we find results consistent within 0.03 dex. 

This consistent evolution in BC03 and PEGASE2 is encour- 
aging. However, to be much more conservative in our esti- 
mate of the uncertainty in the M/L ratio evolution, we take 
the scatter of the models around the best-fit line as our uncer- 
tainty. In Figure 1 1 , in the color range of interest, the scatter is 
approximately ±0.08 dex (20%) at both low and high redshift. 
We use this as the systematic uncertainty in the M/L ratio for 
the purpose of comparing SN rates at low and high redshift in 
§7.3 and §7.4. The uncertainty in the absolute M/L ratio is 
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Figure 11. Evolution of M/L ratio versus color with redshift. Left panel: M/L ratio as a function of u — g color at z = and at z = 1.2 (typical redshift 
in this study). The grid of points show PEGASE2 models with exponentially-decreasing star formation rates with e-folding times r and metallicities Z. For each 
model, star formation begins at z = 4. Models with constant metallicity are connected by solid black lines and models with identical star formation histories are 
connected by dotted lines. For example, models with r = 0, corresponding to a simple stellar population, are the rightmost points (corresponding to Z = 0.01, 
0.02, 0.05) connected by dotted lines. As the models are evolved back in time from an observed redshift of z = to an observed redshift of z = 1.2, the M/L 
ratio decreases and moves away from the Bell03 relation (solid grey line). The dashed grey line shows the relation used in this study for z = 1.2. At z = 1.2 the 
offset from the Bell03 relation is —0.36 dex, or a factor of 0.43. Right panel: Same as left panel, but for g — r color and for an observed redshift of z = 0.6, the 
typical redshift in the rate study of Sharon et al. (2010). The offset here is only —0.14 dex, or a factor of 0.72. 



much greater, due mainly to the uncertainty in the true IMF. 

6. RESULTS AND SYSTEMATIC UNCERTAINTIES 

Here we present our results for the full cluster rate and for 
two galaxy subsets (§6.1) and summarize contributions to the 
uncertainty (§6.2) in each. In §6.3 we show that the rate result 
in the subsets are not sensitive to the specific parameters used 
to select the subset. 

6.1. Results 

The results are presented in Table 7. We derive a rate in 
the full cluster, in red- sequence galaxies only, and in red- 
sequence early-type galaxies only. Each subset includes a 
different number of SNe: As discussed in §3.4, we have 
discovered 8 ± 1 cluster SNe, where the quoted uncertainty 
is due to classification uncertainty (including uncertainty in 
both SN type and cluster membership). Limiting the sam- 
ple to only SNe discovered in galaxies included in the red- 
sequence subset excludes SN SCP06F12 and SN SCP06C1, 
leaving 6.5 ± 0.5 cluster SNe la. The uncertainty here comes 
from the uncertainty in the cluster membership and type of 
SN SCP06E12, which we count 0.5 ± 0.5 cluster SNe la. Fur- 
ther limiting the sample to only SNe discovered in galaxies in- 
cluded in the red-sequence early-type subset, SN SCP06E12 
is eliminated as its host galaxy is dimmer than the ZS50 = 24 
cutoff used for this subset leaving 6 SNe la with negligible 
classification error. The number of SNe la discovered in each 
subset, including classification error, is summarized in Table 7 
under N S n la- 

We normalize the rate in three different ways: by 5-band 
luminosity, by g-band luminosity, and by stellar mass. For 
each cluster, we use the visibility time map T(x,y) (e.g., 
Fig. 7) and the measured luminosity (or mass) profile to carry 
out the integral in equation (2) giving the time-luminosity 



searched. The sum of these values for all 25 clusters is the de- 
nominator of equation (1), the total time-luminosity searched 
in all clusters. This is shown in Table 7 under "Denom" for 
each sample. The rate is simply i a divided by "denom," 
as in equation (1). The contributions to the statistical and sys- 
tematic errors are summarized in Table 8. 

The weighted-average redshift, z, for each subsample is 
given by 

T,iJ x , y T i( x >y) L i( x ,y) 

where Li and Ti are the redshift, luminosity and effective 
visibility time of the i-th cluster, respectively. The weighted- 
average redshift is slightly smaller for the red- sequence and 
red-sequence early-type galaxy subsets. This is because in the 
higher-redshift clusters, a smaller fraction of galaxies meet 
the subset requirements (see z < 1.2 versus z > 1.2 average 
cluster luminosity in Table 6). 

6.2. Summary of Systematic Uncertainties 

Throughout the paper, we have highlighted and addressed 
possible sources of systematic uncertainty. Here we summa- 
rize these sources. In Table 8 we show the relative contri- 
bution of each to the total systematic error, and compare to 
sources of statistical error. 

(1) SN type classification: The uncertainty in the number 
of SNe observed in each galaxy subset was addressed in §6.1. 
The fractional error in the rate is simply the fractional error in 
the number observed. 

(2) Control time: Varying Mb: In our control time simula- 
tions, we assumed a distribution of SN Ia light curve shapes 
and absolute magnitudes. To first order, the impact of these 
assumptions on the control time is captured by varying the 
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Table 7 

Results 



Environment 



Unit 



Nsn la Denom Rate (stat) (sys) 



Full cluster 


SNuB 


1.14 


8.0 ± 1.0 


15.87 


0.50 


+ 0.23 
-0.19 


+ 0.10 
-0.09 


Full cluster 


SNug 






15.96 


0.50 


+ 0.23 
-0.19 


+ 0.10 
-0.09 


Full cluster 


SNuM 






22.41 


0.36 


+0.16 
-0.13 


+0.07 
-0.07 


Red-sequence 


SNuB 


1.13 


6.5 ±0.5 


11.95 


0.54 


+0.25 
-0.19 


+0.07 
-0.07 


Red-sequence 


SNug 






12.20 


0.53 


+0.24 
-0.19 


+0.07 
-0.07 


Red-sequence 


SNuM 






17.61 


0.37 


+0.17 
-0.13 


+0.05 
-0.05 


Red-sequence early-type 


SNuB 


1.10 


6.0 ±0.0 


7.29 


0.82 


+0.39 
-0.30 


+0.09 
-0.08 


Red-sequence early-type 


SNug 






7.59 


0.79 


+0.38 
-0.29 


+0.09 
-0.08 


Red-sequence early-type 


SNuM 






11.77 


0.51 


+0.24 
-0.19 


+0.06 
-0.05 



Note. — "Denom" is 
1O 12 L , B years, 10 12 L©, 9 
and SNuM respectively. 



the denominator of equation (1) and has units of 
years and 10 1 2 Mq years for rate units of SNuB, SNug 



Table 8 

Sources of Uncertainty 



Source of error 


Full 

cluster 

(%) 


Red- 
sequence 

(%) 


Red-sequence 
early-type 

(%) 


Statistical 


Poisson 

Luminosity (stat) 
Luminosity (cosmic var.) 


+40 
-32 

±12 
±16 


+45 
-35 

±6 
±4 


+47 
-36 

±6 
±3 


Total statistical 


+45 
-38 


+46 
-35 


+48 
-37 


Systematic 


SN type classification 
Control time: varying Mb 
Control time: dust distribution 


±13 

+ 8 
-6 

+ 10 

-2 


±8 

+ 8 
-6 


+ 8 
-6 


Luminosity: MAG.AUTO corr. 
Luminosity: K-correction 
Luminosity: Faint galaxy corr. 
Luminosity: r > 0.6(0.8) Mpc 


±7 
±3 

+ 4 
-9 

±4 


±7 
±3 

±1 


±7 
±3 

±1 


Total systematic 


+ 20 
-19 


+ 14 
-12 


+n 

-10 


Total statistical + systematic 


+49 
-42 


+48 
-37 


+49 
-38 



assumed SN la absolute magnitude (§4.3). Variations of 
±0.2 mag resulted in a rate change of tg% 

(3) Control time: dust distribution: In §4.3 we assessed 
the impact of varying amounts of dust extinction on the con- 
trol time. Assuming an unrealistically large amount of dust- 
affected SNe decreased the control time by 9% (increasing 
the SN rate by 10%), while decreasing the amount of dust- 
affected SNe increased the control time by 2% (decreasing 
the SN rate by 2%). We do not apply this systematic error 
to the red-sequence or red-sequence early-type subsets, as we 
have independent evidence that the amount of dust is limited 
in these environments. 

(4) MAGJdJTO correction: In computing the total zsso 
luminosity of each galaxy, we made a correction to the 
MAG_AUTO magnitude ranging from -10% at z 8 50 = 20 
to —30% at ^850 = 25. Varying the range of n used in the 
simulation by ±1 affects the correction by ±7%. 

(5) K -correction: In §5.3, we noted that the scatter of BC03 
templates about the best-fit if -correction is typically less than 
0.03 mag. We use this value as the systematic error on the 
i\~-correction. 



(6) Faint galaxy correction: The average correction C re- 
ported in Table 5 is 1.054. Varying Mg by ± 0.5 magnitudes 
results in an average correction of 1.032 and 1.092 for —0.5 
and +0.5 magnitudes, respectively. Varying a by ±0.2 results 
in an average correction of 1.027 and 1.098 for a = —0.7 and 
— 1.1, respectively. Concurrently varying M% and a within 
the same ranges results in a minimum average correction of 
1.015 (Mb = -22.2, a = -0.7) and a maximum average 
correction of 1.154 (M% = -21.2, a = -1.1). Conser- 
vatively, we assign as the systematic error on the rate 
associated with this correction. This error is not applied to 
the red-sequence or red-sequence early-type subsets because 
these subsets do not include light from galaxies below the de- 
tection threshold. 

(7) Luminosity at large radii: In §5.5 we assumed a model 
for the cluster luminosity profile at r > 0.6 Mpc (0.8 Mpc for 
red-sequence and red-sequence early-type subsets). Varying 
the model luminosity by ±20% resulted in a ±4% change 
in the full cluster rate. The change is much smaller (±1%) 
for the red galaxy subsets because the model is only used at 
r > 0.8 Mpc. 

(8) M/L ratio: In §5.7 we used a M/L ratio to translate 
stellar luminosity to stellar mass. Rather than estimating the 
absolute uncertainty in the M/L ratio (which is strongly de- 
pendent on assumptions), we estimated the uncertainty in the 
evolution of the M/L ratio from low to high redshift. This is 
the relevant uncertainty for comparing rates at differerent red- 
shifts in order to derive the SN la delay time distribution. We 
defer discussion of this uncertainty to §7.4 where we discuss 
uncertainties in the DTD. 

6.3. Effect of Varying Subset Requirements 

In selecting our red- sequence and red- sequence early-type 
galaxy subsamples, we required red-sequence galaxies to be 
within ±0.2 mag of the color of their cluster red sequence. 
For early-type galaxies, we required the asymmetry parame- 
ter to be < 0.1 and the Gini coefficient to be > 0.40. It is 
interesting to test the sensitivity of the results to variations in 
the requirements. In Figures 12 and 13 we vary the require- 
ments and observe the effect on the rates. As requirements are 
made more strict (for example, narrowing the red sequence) 
the total mass of the sample decreases. At the same time, 
SNe fall out of the sample when their host galaxies are cut. 
The Poisson error increases as the number of included SNe 
shrinks. 



The SN la Rate in High-Redshift Galaxy Clusters 



25 




0.0 

0.05 0.10 0.15 0.20 0.25 
Red sequence half-width (mag) 

Figure 12. The effect of varying the width of the red sequence. The nominal 
red-sequence rate result corresponds to a half-width of 0.20 mag. The inner 
and outer error bars represent the statistical and total uncertainty, respectively. 
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Figure 13. The effect of varying the morphology parameter requirements. 
Negative A values correspond to a more strict selection and a higher-purity 
early-type galaxy sample. The requirements are asymmetry < 0.1 + A and 
Gini coefficient > 0.40 — A. The nominal red-sequence early-type rate 
corresponds to A = 0. The red-sequence half-width is fixed at 0.2 mag. 
The inner and outer error bars represent the statistical and total uncertainty, 
respectively. 

There is not a strong dependence of the SN la rate with 
galaxy color residual from the red sequence (Fig. 12). Even 
in cluster galaxies that lie in a tight range around the red- 
sequence (±0.08 mag), we find a SN la rate consistent with 
the full cluster rate. Similarly, there is no significant rate trend 
with the purity of the early-type sample (Fig. 13). We hap- 
pened to pick morphology requirements that yield a slightly 
higher rate than other choices, but such variations are ex- 
pected with small-number statistics and are accounted for by 
the Poisson uncertainty in the result (Tables 7 and 8). Even in 
the most-selective subset (A = —0.04), the rate is consistent 
with the full cluster rate. 

7. DISCUSSION AND CONCLUSIONS 

7.1. Host-less Cluster SNe la 

As reported by Dawson09, we have discovered one poten- 
tial host-less cluster SN la among the 8 ± 1 cluster SNe la. 
SN SCP06C1 is projected near two possible host galaxies: A 
zs5o = 21.6 spiral galaxy I'M West of the SN, and a sig- 
nificantly fainter zs5o = 24.6 galaxy 0".45 (~3.5 kpc at the 



cluster redshift) Northeast of the SN (See Dawson09, Fig. 2). 

The galaxy- subtracted SN spectrum clearly shows a SN la 
at redshift z = 0.98 near maximum light, consistent with the 
light curve fit. The redshift ofz = 0.98±0.01is consistent 
with the cluster redshift of 0.974. The bright spiral galaxy is 
actually in the background of the cluster, at z — 1.091. Strong 
[Oil] emission is visible in the spectrum, along with Ca H & 
K and US absorption. Unfortunately, the small separation be- 
tween the main galaxy and the smaller galaxy to the Northeast 
means that the spectrum of the smaller galaxy is dominated by 
light from the larger galaxy, making it impossible to assess a 
redshift. It is thus possible that the small galaxy is at the clus- 
ter redshift and is the actual host of the SN. Alternatively, the 
small galaxy might be at the same redshift as the larger galaxy 
and physically associated with it (either as a satellite galaxy 
or as part of the spiral structure of the galaxy). It is interest- 
ing to note that the SN is only 20" (160 kpc) projected radius 
from the center of the cluster, perhaps giving more weight to 
the hypothesis that it is associated with a diffuse intracluster 
stellar component. 

Not being able to confirm or reject this SN as host-less, 
we have an upper limit of one host-less SN out of a total of 
8 ± 1. Discovering one host-less SNe la out of seven total 
would imply an intrinsic host-less SN la fraction of 14%+^T° 
(binomial 68% confidence intervals), and a 95% upper limit of 
< 47%. This is broadly consistent with host-less SN la con- 
straints at intermediate redshifts (Sharon et al. 2010) and at 
low redshift (Gal- Yam et al. 2003; Sand et al. 2011). At low 
redshift it has been possible to confirm the host-less nature of 
some SNe using deeper follow-up imaging, leading to better 
constraints. The upper limit of < 47% is also consistent with 
direct measurements of intracluster light at low redshift, but 
does not strongly constrain evolution. A sample twice the size 
or larger, with deeper follow-up to confirm host-less SNe la 
would begin to place interesting constraints on hypotheses for 
the formation of the intracluster stellar component from z > 1 
to today. 

7.2. Comparison to Other Cluster Rate Measurements 

Cluster SN la rates have been reported at lower redshifts 
by several groups. In nearby (z < 0.2) clusters, measure- 
ments include those of Sharon et al. (2007) at z ~ 0.14, 
Mannucci et al. (2008) at z ~ 0.02, and Dilday et al. (2010) 
at z ~ 0.09 and z ~ 0.22. At intermediate redshifts, 
Sharon et al. (2010) recently reported the rate in 0.5 < 
z < 0.9 clusters (median z ~ 0.6). At higher redshifts, 
Gal- Yam et al. (2002) placed the first constraints on the z > 
0.8 cluster rate using a sample of three clusters at z = 0.83, 
0.89 and z = 1.27. However, their SN sample included only 
one firm SN la at z = 0.83. The resulting rate has corre- 
spondingly large uncertainties and essentially places only an 
upper limit on the z > 0.9 cluster rate. Our result is thus a 
large step forward in the measurement of the SN rate in the 
highest-redshift clusters. 

In Figure 14 we compare our full cluster rate to the 
lower-redshift rate measurements that have been normalized 
by stellar mass, permitting a comparison across redshifts. 
Here we have made an adjustment to the value reported by 
Sharon et al. (2010). Sharon et al. used the mass-to-light ra- 
tio of Bell03 for the SDSS g and r bands, but did not apply 
a correction for evolution between z ~ 0.6 and z = 0. Us- 
ing the method described in §5.7 we find that a —0.14 dex 
offset should be applied to the mass to account for evolution 
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Figure 14. Cluster rate measurements (all galaxy types) from this work and 
the literature. The rate of Sharon et al. (2010) shown has been adjusted up- 
ward by 38% from the reported rate (see text). The top axis shows the 
time elapsed since an assumed cluster formation redshift of Zf = 3. The 
solid grey line represents the SN la rate for the best-fit power-law DTD: 
iaM — ^(t)/m(t), where ^(t) oc t s . The dotted grey lines show 
the range of la error on s. 



from z = 0.6 to z = (Fig. 11, right panel). We therefore 
adjust the reported rate of Sharon et al. upward by 0.14 dex 
(38%). The rate compilation of MaozlO reflects this adjust- 
ment. Whereas the adjusted Sharon et al. rate shows an in- 
dication that the cluster rate is increasing with redshift, for 
the first time we find an increasing rate with high significance 
(> 2a). 

We point out that the popular "A + £>" model 
(Scannapieco & Bildsten 2005) is insufficient for describing 
the change in cluster rate with redshift. In the A + B model 
the SN rate is the sum of a term proportional to the total stel- 
lar mass and a term proportional to the recent star formation 
rate: 7^sn ia = AM* + BM*. This simple model is conve- 
nient for predicting the SN rate in environments with vary- 
ing amounts of recent star formation as it accounts for the 
increased SN Ia rate at short delay times. (In fact, we use 
this model in Meyers 1 1 to derive limits on the expected ratio 
of SNe Ia to SNe CC in early-type galaxies.) However, the 
model lacks theoretical motivation and breaks down in other 
situations. For example, Greggio et al. (2008) note that it can- 
not adequately describe the observed contribution from SNe 
with intermediate delay times (e.g., Totani et al. 2008). This 
point is reinforced by the observation of a changing cluster 
rate with redshift: In clusters, the A component is dominant 
at all redshifts observed. As M* is not changing significantly 
with redshift, the rate would be expected to remain constant 
under this model. Instead, we require a DTD model wherein 
the rate decreases at large delay times (as it does in most the- 
oretically motivated models). 



7.3. The Cluster SN Ia Delay Time Distribution 

The cluster rates constrain the SN Ia delay time distribu- 
tion, ^(t), over the range of delay times from a few Gyr to 
~ 10 Gyr. To illustrate the cluster rate constraints, we pa- 
rameterize the DTD with a power law in time: ^(t) oc t s . 
A power law is not only a convenient parameterization in the 
face of limited data, but is a theoretically motivated function 
for the DD scenario, where the late- time (t > 1 Gyr) DTD 
shape is set by the distribution of WD separation after the sec- 
ond CE phase and the merger timescale due to gravitational 
radiation (Greggio 2005). 

We make the approximation that all clusters formed in a 
single burst of star formation at Zf = 3 and that the age 
of the stellar population therefore corresponds to the elapsed 
time from Zf to the cluster redshift (Fig. 14, top axis). While 
clearly a simplification, a single star-formation burst captures 
the idea that the timescale over which star formation oc- 
cured in cluster early-type galaxies is short compared to the 
time since star formation ceased. The assumed burst redshift 
Zf = 3 is consistent with measurements of cluster early-type 
galaxies showing that star formation was mostly completed by 
this redshift (e.g., Gobat et al. 2008). Below, we show that the 
derived DTD is relatively insensitive to the redshift assumed. 

As noted in §5.7, the DTD is normalized by initial stel- 
lar mass, whereas the cluster rate measurements (including 
ours, for consistency) have been normalized by current stel- 
lar mass. The DTD, Sfr(t), is therefore related to the cluster 
rate by &(t) = m(t)7ZsN i&(t) where m(t) is the fraction of 
stellar mass remaining at time t after the star formation burst. 
The specific choice of m(t) does not have a significant impact 
on the derived DTD: regardless of the model or IMF assumed, 
the stellar mass declines by only ~10% over the age range of 
interest, ~ 3 to 11 Gyr. For consistency with MaozlO, we 
use the remaining stellar mass fraction tabulated by BC03, 
^bco3(^), but corrected to m(t) = 1 — (1 — mBC03(^))/0.7 
to effectively convert from the Salpeter IMF used in BC03 to 
a "diet" Salpeter IMF. This correction has only a very small 
effect on the result (see below). 

We find a best-fit value of 

s = -1.41™ (13) 

using the statistical+ systematic error (added in quadrature) 
reported for each rate measurement. In Figure 14, the 
solid grey line shows the best-fit cluster rate for this value: 
T^SNiaO) = #(t)/ra(t), where ex t~ 1A1 . Note that 
the x 2 of me best-fit model is surprisingly small: 0.40 for 4 
degrees of freedom. The a priori probability of finding a \ 2 
smaller than 0.40 is less than 2%. This is difficult to under- 
stand given that the measurement errors are generally domi- 
nated by Poisson noise in the number of SNe observed and 
are thus unlikely to be overestimated. 

The best-fit value is consistent with measurements of the 
late-time DTD in the field (Totani et al. 2008). Most pre- 
dictions for the SD scenario show a steeper late-time DTD 
(Greggio 2005; Ruiteretal. 2009; Mennekens et al. 2010) 
with an effective value for s ranging from s ~ —1.6 (Greggio 
2005) to s < —3 (Mennekens et al. 2010), depending on the 
details of the scenario and binary evolution. However, some 
groups have found that the SD scenario could be consistent 
with a less-steep DTD (s ~ — 1) given the right combination 
of main sequence and red giant secondaries (Hachisu et al. 
2008). In the DD scenario, the predicted shape of the DTD de- 
pends on the distribution of binary separations after the com- 
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mon envelope phase of the WDs, a difficult distribution to 
predict. However, a slope of s = —1.4 (and a range of similar 
values) would not be surprising in the DD scenario. 



7.4. Additional DTD Systematic Uncertainties 

Variations in the assumed cluster star formation, initial 
mass normalization and mass-to-light ratio evolution have a 
small affect on s compared to the measurement error. 

(1) Age of clusters' stellar populations: Above, we as- 
sumed a single burst of star formation at Zf = 3. Moving 
this single burst to Zf = 4 results ins = —1.55. A more 
recent burst, Zf = 2.5, results in s = —1.30. MaozlO give 
a treatment of variations from the single-burst approximation, 
also finding that the affect on s is small. 

Our rate measurements in red and early-type galaxies pro- 
vide a good consistency check that recent star formation does 
not significantly contribute to the SN la rate: if it did, we 
would observe a higher rate in the full cluster than in these 
subsamples. Surprisingly, we observe the opposite trend (al- 
though the significance is low). The red- sequence early-type 
subsample includes 53% of the stellar mass of the full clus- 
ter sample, and 6 SNe la. The remaining 47% of the full 
cluster sample (which includes bluer galaxies and late-type 
red- sequence galaxies) accounts for only 2 ± 1 SNe la. At 
low redshift, Mannucci et al. (2008) found a similar trend be- 
tween E/SO galaxies and SOa/b galaxies within 0.5 Mpc of 
cluster centers, though also at < la significance. 

(2) Remaining stellar mass: Whereas the DTD is normal- 
ized by initial stellar mass and cluster rate measurements have 
been normalized by current stellar mass, we have assumed 
a remaining stellar mass fraction m(t) to convert from cur- 
rent to initial stellar mass. Although different models and 
IMFs can yield sigificantly different m(t), we are only con- 
cerned here with the change in m(t) between ~ 3 Gyr and at 
~ 11 Gyr. (The absolute value of m(t) affects only the nor- 
malization of ^(t), with which we are not concerned.) For- 
tunately, the evolution in m(t) in this age range is small and 
consistent between models, and so the effect on s is small. For 
example, using mBC03 (t) (assuming a Salpeter IMF) rather 
than correcting to a diet Salpeter IMF (as we have done) only 
changes the best-fit value from s = — 1.41 to s = —1.38. 

If in §5.7 we had used a M/L ratio directly normalized by 
initial mass, rather than normalizing by current mass and later 
converting to initial mass, the results would be very similar. 
(We have not done this for consistency with other rate mea- 
surements.) In the PEGASE2 models in Figure 11 (left panel) 
evaluted at z = 1.2, the ratio of current to formed stellar mass 
varies slightly across the models, but is fully contained in the 
range 0.66 ± 0.03. The same models evaluted at z = have 
a ratio of 0.59 ± 0.03. This is consistent with the ~ 10% 
evolution in m(t) over this range as tabluated by BC03. 

(3) M/L ratio evolution: While the overall normalization 
of the M/L ratio will only affect the normalization of ^(t) 
and not s, the evolution of the M/L ratio will affect s. In 
§5.7 we assigned a liberal 20% systematic uncertainty to the 
evolution of the M/L ratio over the redshift range of interest. 
To estimate the effect of this systematic uncertainty, we adjust 
our rate measurement by 20% and that of Sharon et al. (2010) 
by 10% and refit s. The resulting change in s for positive and 
negative shifts is —0.15 and +0.18 respectively, less than half 
of the nominal error in s. 



7.5. Conclusions 

In this paper, we have made a measurement of the high- 
redshift cluster SN la rate. Thanks to an unusually com- 
plete dataset (particularly for a rate study) the measurement 
is quite robust, with statistical and systematic uncertainties on 
par with or better than measurement uncertainties at low red- 
shift. We highlight several important and/or unique aspects of 
the measurement: 

• The SN classification approach takes advantage of all 
relevant information. Thanks to the "rolling search" 
strategy of the survey and the nearly complete spec- 
troscopic follow-up, most candidates have a full light 
curve and a host galaxy redshift, greatly reducing clas- 
sification uncertainty. 

• The position-dependent control time allows one to cal- 
culate a supernova rate given an arbitrary observing pat- 
tern and luminosity distribution. 

• The control time calculation includes a full distribution 
of SN properties and the systematic uncertainty associ- 
ated with the assumed distribution is carefully quanti- 
fied. Thanks to the depth of the observations, the de- 
tection efficiency approaches 100% during the period 
of the survey for most of the clusters, meaning that the 
systematic uncertainty is low. 

• Statistical uncertainties associated with the cluster lu- 
minosities, including both statistical variations and cos- 
mic variance, are included in the total uncertainty. 
Also, light in the outskirts of each galaxy (outside the 
S Extractor MAG_AUTO aperture) is accounted for. 
This is a significant component of the total cluster lu- 
minosity. 

• Cluster SN la rate measurements are normalized con- 
sistently across redshifts using a redshift-dependent 
mass-to-light versus color relation. 

For the first time our result shows at the > 2a level that the 
cluster SN la rate is increasing with redshift. Simply by com- 
paring the low- and high-redshift cluster rate measurements, 
the shape of the late-time SN la delay time distribution can be 
constrained. The power of the measurement for this purpose 
comes both from the high redshift and relatively low statis- 
tical and systematic uncertainties in the measurement. While 
we cannot conclusively rule out either the single degenerate or 
double degenerate class of progenitors via the delay time dis- 
tribution, the binary evolution that could lead to each model 
are constrained. The DD scenario is consistent with the mea- 
surement under a wide range of plausible binary evolution pa- 
rameters, while there is a stronger constraint on binary scenar- 
ios that could lead to an SD scenario. Finally, this measure- 
ment is unique in constraining the delay time distribution at 
delay times of a few Gyr. In future studies, it can be used in 
combination with other cluster rates and other delay time dis- 
tribution measurements (e.g., MaozlO) to place even tighter 
constraints on models for binary evolution and SN la progen- 
itor scenarios. 
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